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METHOD AND APPARATUS FOR READING 
REPORTER LABELED BEADS 
Related Applications 

\ *j This application is based on prior copending provisional application Serial 

rf| 5 No. 60/240,125, filed on October 12, 2000, and prior copending provisional application 

111 Serial No. 60/242,734, filed on October 23, 2000, the benefits of the filing dates of 

I -j which are hereby claimed under 35 tlS.C. § 1 19(e). This application is further based 

.. * on prior copending conventional application No, 09/939,292, filed on August 24, 2001 , 

M which is based on a prior copending provisional application Serial No. 60/228,076, 

; s, f 10 filed on August 25, 2000, the benefits of the filing dates of which are hereby claimed 

nx under 35 U.S.C. § 1 19(e). and under 35 U.S.C. § 12a 

Field of the Invention 
The present invention generally relates to a method and apparatus employed to 
image an encoded reporter labeled bead, and more specifically, to a method and 
15 apparatus that enable encoded reporter labeled beads to be imaged in stasis or when 
entrained in a flow of fluid, such that imaging data thus obtained can be employed to 
decode each encoded bead by determining the identity of reporters bound to that bead* 

Background of the Invention 
Single nucleotide polymorphisms (SNPs) are locations in the genome where a 
20 single base substitution has occurred. SNPs are estimated to occur in the human 
genome at a frequency of approximately 1:1000, implying that there are three million 
SNPs in the three billion nucleotide human genome. Since most known gene 
sequences are on the order of 1,000 base pairs in length, each gene is expected to 
contain one SNP. It is estimated that there are 100,000 human genes, meaning that 
25 there are 100,000 SNPs, which may directly affect the function and/or expression of the 
resulting proteins. 

The relative abundance of SNPs in the genome has stimulated efforts to 
quantify the location and frequency of occurrence of single base substitutions as a tool 
for the analysis of gene function. Methods for the detection of SNPs include the 
30 oligonucleotide ligation assay (OLA), single-strand conformation polymorphism 
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analysis, allele-specific oligonucleotide (ASO) hybridization, and the single base chain 
extension (SBCE) assay. 

While SNPs are abundant, an individual SNP is not rich in information. Most 
SNPs are in non-coding or non-regulatory regions of the genome and may not affect 
5 gene expression at all Of those SNPs that occur within the coding regions, many are 
likely to occur in non-critical regions of the resulting protein or may result in a benign 
(or nonexistent) amino acid substitution and therefore have little or no bearing on the 
protein's function. Further, only a fraction of the total number of genes are actively 
expressed at a given time, so the presence of an SNP within a gene does not indicate a 
10 priori that it has significant phenotypic relevance. Since the fraction of phenotypically 
iJ3 relevant SNPs is small, it is useful to have a high throughput means of analyzing SNPs 

; 3 in order to identify those of biological importance. Once the genome has been 

I jj thoroughly analyzed and the locations and relative abundances of important SNPs have 

ry been identified, a multiplexed method of SNP analysis from an individual's 

1 5 deoxyribonucleic acid (DNA) will be clinically useful. It would therefore be desirable 
to have a high throughput method for analyzing numerous SNPs in a short period of 
* time. 

The locations of SNP and other polymorphic loci within the genome can be 
determined by sequencing and comparing the same sections of the genome from 
20 numerous individuals. Those locations within the genome that are statistically variable 
across individuals are polymorphic. The biological relevance of a given polymorphism 
can be determined by correlating the different alleles to the presence of disease or other 
phenotypic traits. Hence, there is a need for a robust, inexpensive, and widely available 
method for sequencing gene-sized lengths of DNA in order to discover the locations 
25 and biological relevance of polymorphic sites. 

One SNP analysis method calls for the binding of oligonucleotides to supports 
such that numerous identical oligos are bound to a solid support, and so that different 
supports bear different oligo sequences. One method of encoding an oligonucleotide 
library useful for SNP analysis is to place unique optical reporters on solid supports 
30 during combinatorial chemical synthesis. The attachment of reporters to the supports 
may be by means of covalent bonds, colloidal forces or other such means, to ensure that 
reporters stay in contact with, or in close proximity to, the solid support The solid 
supports are typically beads of polystyrene, silica, resin, or any another substance on 
which compounds can be readily synthesized and to which reporters can be affixed in a 
35 spUt/add/pool (SAP) combinatorial process. Each reporter encodes both the identity of 
a molecular component as well as its place in the synthetic process. By enumerating 
the optical characteristics of each reporter bound to a solid substrate it is possible to 
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decode libraries of unique compounds numbering in the billions. As noted above, 
useful genetic assays can be performed by combinatorially synthesizing 
oligonucleotides on a bead Hbraiy such that a given bead bears numerous identical 
covalently bound oligos and each bead in the library bears a different oligo sequence. 
5 In addition to its oligo sequence, each bead bears a unique optical signature comprising 
a predefined number of unique reporters, where each reporter has a predefined 
combination of different fluorochromes. A bead's optical signature is correlated to the 
addition sequence of each reporter during the synthetic process to enable identifying the 
unique nucleotide sequence on that bead. By imaging the beads, the optical signatures 
10 can be read and correlated to the corresponding oligo sequences. 
% In addition to DNA analyses, the reporter labeling method is also useful for 

synthesizing diverse libraries of chemical compounds on beads for subsequent analyses 
in dnig candidate screening. Likewise the same method can be used for synthesis of 
protein libraries on beads where the base unit of synthesis is one of twenty amino acid 
15 sequences. 

Generally, with only a few reporters and colors, the number of unique 
signatures that can be created is quite substantial. For example, using only five colors 
and five reporters, more than 10,000 unique signatures can be generated. Using six 
colors and 10 reporters, over 115 million unique signatures can be generated to create a 
20 very diverse bead library. Clearly, the number of unique combinations that can be 
identified using reporter labeled beads is substantial. The nature of the apparatus 
required to identify the unique spectral signatures of such beads is discussed in greater 
detail below. 

Substrate-Based Approach to Analyzing Reporter labeled Beads 
25 In order to read the reporter signature of a bead, an image of the bead must be 

acquired with sufficient spatial resolution to discriminate the locations of individual 
reporters. If reporter size or shape are used in the signature scheme, the spatial 
resolution must be sufficient to &scniiiinate these parameters as well. Further, the 
acquired image must have sufficient spectral resolution to accurately discriminate the 
30 multiple colors emitted from a single reporter. Further still, the quality of the imagery 
of each bead acquired must be sufficient to ensure that at least one copy of every unique 
reporter on a bead is evident in the view. Even when multiples of each reporter are 
bound to a bead, there remains a probability that not every unique reporter will be 
resolved in a given image. Reporters may not be in clear focus or they may not be 
35 exposed to the optical collection system because of their disposition on the bead. In 
such cases, multiple images should be acquired of each bead at different focal planes or 
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from different perspectives to ensure at least one copy of every unique reporter is 
successfully discriminated 

One technique that might be used to read the reporter signature on beads using 
a conventional fluorescence microscopy apparatus would require that the beads be laid 
down on a planar substrate in order to present an optically readable bead array. If the 
beads were used in an assay prior to being affixed to a substrate, the assaying process 
should not disrupt the bound reporter signature. In the case where beads are affixed to a 
substrate before being used in an assay, the advantageous kinetics associated with the 
large surface area to volume ratio of free beads is lost Nevertheless, the limitations of 
conventional microscopy techniques impose the requirement that beads be affixed to a 
substrate prior to analysis, thereby adding numerous preparatory steps to bead-based 
assays. These preparatory steps add time and expense, while simultaneously reducing 
the flexibility and utility of bead-based analytical processes. 

In one bead arraying method, described in U.S. Patent No. 5,855,753, beads are 
placed on a substrate and caused to form a fixed monolayer through the use of an 
electric field. An "electrochemical sandwich" is formed by suspending the beads in an 
electrolytic fluid placed between an anode and a cathode. Using either an alternating 
current (AQ and or a direct current (DC) field applied to the sandwich in an 
appropriate manner, over time, the beads are caused to aggregate in specific groups, in 
a monolayer on the substrate. 

Another bead arraying method is described in International Patent 
Application WO97/40385, which indicates that the electrochemical sandwich method 
is further enhanced by use of a specialized electrode in conjunction with externally 
applied illumination patterns that serve to further control the electrokinetic forces, 
which mediate bead aggregation on a substrate. U.S. Patent No. 5,695,934 discloses 
yet another method in which beads are laid on a substrate and affixed by chemical 
affinity between the functionalized surface of the substrate and bead-bound moieties. 
Other methods for arraying beads on substrates exist, but in all such methods, the goal 
is to ensure that the bead layer is fixed in place, preventing movement of the beads 
during the process of reading the beads. Typically, in any microscopy process for 
reading beads, it would be necessary to affix the bead-substrate to a two-axis stage and 
move the stage in a pattern that enables each portion of the substrate to be read. This 
process involves numerous cycles of acceleration and deceleration of the substrate as it 
is moved on the stage, which would likely induce independent movement of the beads, 
if they are not securely affixed to the substrate. 

Bead movement is not the only complication associated with reading labeled 
beads on a planar substrate. Another consideration is the need for achieving an 
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accurate focus across the field of view (FOV), which can be compromised by any 
non-planarity of the packed beads on the substrate or by any non-planarity of the 
substrate itself. For these reasons, the focus on each portion of the bead array must be 
individually achieved to ensure proper resolutioa Although autofocus systems are well 
known in the art, the focus step requires additional time, expense, and adds variability 
to the process. Additional images may be required to discriminate the different 
fluorescence emission spectra of bead-bound signaling molecules and of the reporters 
themselves. In addition, if the signaling molecules or reporters are randomly 
distributed on the beads, it may not be possible to identify signals or signatures from a 
fraction of the beads due to the absence of the signals or reporters from the planar FOV. 
The planar substrate preparation presents only one, or at most two, of six possible 
perspectives from which to view the bead, increasing the likelihood that reporters will 
be hidden from the imaging system. 

The complexity associated with arraying beads hinders bead-based analytical 
approaches, regardless of the number of beads in the library. As the size of an 
analytical bead library grows beyond roughly a million beads, the substrate-based 
approach to bead imaging becomes highly impractical. Significant difficulties arise 
when tens of millions to billions of beads must be analyzed, requiring that the bead 
array substrates grow substantially in size. Hie difficulties involved in creating a 
uniform, tightly packed, fixed array increase greatly with the size of the array. 
Furthermore, accurate and rapid positioning of the array during the imaging process 
becomes far more difficult. The size, expense, and low throughput of such systems rule 
out their widespread use in research and for point-of-care applications. Therefore, an 
improved method for analyzing beads is desired Preferably this method should 
eliminate the need for placing the beads on a substrate and enable the simultaneous 
imaging of multiple focal planes and multiple bead orientations. Ideally, this new 
method would simultaneously image the entire spectrum of bead fluorescent emissions 
and provide enough spectral resolution to discriminate the colors originating from each 
reporter. Finally, an ideal method would conveniendy handle billions of beads and 
enable ultra-high speed imaging to analyze large bead libraries in a matter of hours. 

Summary of the Invention 

In accord with the present invention, a method and apparatus for imaging and 
reading reporter labeled beads are provided. The method and apparatus of the present 
invention enable individual encoded beads to be imaged, and the compound attached to 
that bead to be identified as a function of image data, so that the identity and sequence 
of all sub units of the compound are determined. Generally, a single encoded bead will 
include redundant copies of the same compound and the associated reporters. In one 
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preferred application of the present invention, the sub units are oligos, the resulting 
compound is an oligonucleotide, and each reporter is preferably discriminable by its 
color. However, other types of sub units forming other types of compounds and 
reporters that are uniquely discriminable by characteristics other than color can be 
5 imaged and decoded using the method and system of the present invention. Essentially, 
any characteristic that is determinable via imaging, such as size and shape, can be 
employed as a reporter characteristic. 

It should also be understood that the method and apparatus of the present 
invention may be applied to image and decode reporter labeled beads independently of 
Q 10 the process employed to produce such reporter labeled beads. It is anticipated that 

^ processes such as stochastic synthesis, directed synthesis, combinations thereof, or the 

\* attachment of a pre-synthesized compounds to previously (or subsequently) encoded 

© reporter labeled beads will be employed to produce reporter labeled beads. Unlike a 

! ^ stochastic SAP process, in a directed synthesis the path of each bead is predetermined 

%j 15 during compound synthesis and reporter addition. When reporter labeled beads are used 
« with pre-synthesized compounds all unique reporter types may be bound to the carrier 

in one step to create a unique signature for the bead In this case the unique combination 
U of reporters bound to a bead simply create a unique identity and do not necessarily 

rU encode the subunit sequence for the compound bound to the bead. A cross reference 

£f 20 table or other means may be created to correlate bead signature to compound identity. 

In this case the present invention may be used to identify the unique bead signature to 
determine compound identity by searching the cross reference table. 

The basic method involves focussing light from an encoded bead along a 
collection path, and dispersing the light into a plurality of light beams as a function of a 
25 plurality of different discriminable characteristics of the light that are indicative of the 
plurality of different reporters associated with the encoded bead. Each light beam is 
focussed to produce a respective image corresponding to that light beam. A plurality of 
signals are generated in response to the respective images thus produced. Each signal 
generated indicates whether a different one of the plurality of reporters is associated 
30 with the encoded bead. The signals are analyzed to decode a sequence and to identify 
the sub units forming a component attached to the encoded bead, and these steps are 
repeated for successive encoded beads. 

Preferably, die step of dispersing comprises the steps of dividing the light into 
the plurality of light beams as a function of the wavelength of the light beams. Also 
35 preferably, the step of analyzing the plurality of signals comprises the step of 
constructing a sequence library of the plurality of components based on each encoded 
bead that is decoded. 
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It may be necessary to illuminate an encoded bead with light to enable the 
encoded bead to be read. In one embodiment, the encoded beads are entrained in a 
fluid, and the encoded beads are imaged in a flow of the fluid However, the present 
invention does not require that the beads be imaged in a fluid flow. 

Preferably the encoded beads are analyzed by determining one or more portions 
of each image corresponding to a unique reporter, and then detenmning a signature of 
each reporter based on the image data. The reporters thus identified determine the 
identity and the sequence in which the sub units that form the compound attached to the 
encoded bead were added 

A predefined number of unique reporters should be detected on an encoded 
bead The method can also include the step of disregarding all signals relating to an 
encoded bead if the plurality of signals for that encoded bead indicate that fewer 
reporters are associated with the encoded bead than were expected This step can be 
facilitated by referring to an encoded bead legend that relates each unique reporter to a 
specific sub unit and a specific disposition of that sub unit in a sequence that forms a 
compound 

It should also be noted that encoded beads are usable to analyze DNA 
sequences, or for other DNA related research. In such an application, it is preferable 
not to collect data for encoded beads that have not experienced a binding event Beads 
that have experienced a binding event are referred to herein as ••positive" beads. In at 
least one embodiment of the present invention, only data from positive beads are 
analyzed 

Several different embodiments of imaging apparatus are usable to image and 
decode encoded beads. In general, the apparatus will include a collection lens to collect 
and focus light from an encoded bead in a desired direction, a dispersing component 
that receives the light from the collection lens and disperses the light into a plurality of 
light beams as a function of a plurality different discriminable characteristics of the 
light, such that the plurality of different discriminable characteristics are indicative of a 
plurality of different reporters. Also included is at least one pixelated detector, and an 
imaging lens that focusses each of the plurality of light beams on the at least one 
pixelated detector, thereby producing a respective image corresponding to each of the 
plurality of light beam. The pixelated detector provides an output signal for each 
respective image, such that each output signal may individually, or in combination with 
other output signals, indicate whether a different reporter is present on the encoded 
bead A signal processor is coupled to the pixelated detector to receive its output signal 
and processes the output signal to decode the identity and sequence of the sub units that 
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form the compound attached to an encoded bead, or alternatively, to decode the bead to 
identify the compound using a cross reference table. 

Preferably, in one exemplary application, the signal processor is adapted to 
generate sequence contigs from a plurality of decoded sequences representing a 
plurality of encoded beads. Such sequence contigs preferably are used to identify at 
least one of a genomic DNA sequence, a polymorphic allele, and an expressed gene. 

The signal processor can be adapted to disregard all output signals relating to a 
reporter if signals from an identical reporter have already been analyzed The signal 
processor can similarly be adapted to disregard all output signals for an encoded bead if 
the signals indicate that few than a predetermined number of reporters are associated 
with the encoded bead The signal processor is preferably adapted to employ an 
encoded bead legend that relates each unique reporter to a specific sub unit and a 
specific disposition in a sequence of sub units that comprise the compound associated 
with the encoded bead. By employing such an encoded bead library, the signal 
processor can selectively disregard all output signals if it is determined that the encoded 
bead does not correspond to that encoded bead legend Alternatively, the signal 
processor can be adapted to employ a cross reference table (such as the reporter legend 
of FIGURE 19) to relate bead identity to compound identity. Also, the signal processor 
can be controlled to only analyze signals relating to positive encoded beads. 

In at least one embodiment, the encoded beads are entrained in a fluid and 
imaged while the fluid is flowing through the imaging system. The imaging system 
includes a bead reservoir containing encoded beads that have not been imaged, and 
fluid channels to control a flow of fluid A light source is optionally included to 
illuminate the encoded bead to be imaged. 

The dispersing component is preferably either a prism or a plurality of dichroic 
filters. If a prism is employed, the images will be convolved and the signal processor 
will need to deconvolve the images. One or more time delay and integration (TDI) 
detectors are also preferably included in the imaging system to detect the light from the 
beads. In one embodiment, a plurality of dichroic filters, imaging lenses, and detectors 
are included 

Still another aspect of the present invention relates to a method for 
simultaneously imaging a plurality of reporters disposed on substantially different 
portions of an encoded bead, to identify each unique reporter included on the encoded 
bead The method includes the steps of receiving light from the encoded bead along a 
plurality of collection paths that are substantially spaced apart, such that light from the 
reporters disposed on the different portions of the encoded bead affect the light received 
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thereftom. The light is processed to identify each unique reporter included on the 
encoded bead 

A final aspect of the present invention is directed to a method for employing an 
oligo library encoded on beads for at least one of DN A sequencing, a polymorphism 
5 analysis, and an expression analysis. An imaging system capable of decoding a 
sequence of encoded beads conveyed in a flow of fluid is provided Essentially, a 
complete encoded bead library of N-mer oligos is generated, and the library is used to 
perform either the DNA sequencing, the polymorphism analysis, or the expression 
analysis. Positive beads are imaged and decoded using imaging data produced by an 
10 imaging system such as described above. 

When DNA sequencing is the analysis function that is selected, the genomic 
DNA component is amplified using primers for extended sequences of interest When 
a polymorphism analysis function is selected, the genomic DNA component is 
amplified using primers for polymorphic regions of interest Alternatively, when an 
15 expression analysis function is selected, the ribonucleic acid (RNA) component is 
amplified using primers for genes of interest 

Regardless of the specific analysis function that is selected, the amplified 
component is hybridized to the encoded bead library; and the imaging system is 
employed to identify oligo sequences of encoded beads hybridized as noted above. 
20 From the imaging data, sequence contigs are constructed from the oligo sequences 
identified by the imaging analysis to identify either a genomic DNA sequence, a 
polymorphic allele, or an expressed gene. In a preferred embodiment, the N-mer oligos 
comprise oligos having a length equal to ten. 

Brief Description of the Drawing Figures 
25 The foregoing aspects and many of the attendant advantages of this invention 

will become more readily appreciated as the same becomes better understood by 
reference to the following detailed description, when taken in conjunction with the 
accompanying drawings, wherein: 

FIGURE 1 is an isometric view of an embodiment of an imaging system useful 
30 for imaging reporter beads, in which multiple legs are employed for the spectral 
decomposition and imaging to collect light signals from multiple perspective positions; 

FIGURES 2A and 2C are isometric views of two related embodiments of an 
imaging system using a prism for spectral dispersion of light and usable for imaging 
reporter beads, wherein FIGURE 2C also includes a slit for spatial filtering of 
35 extraneous light; 

FIGURE 2B is an schematic view of embodiment 2A showing bead traversal 
along a direction of motion in both object and image space; 
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FIGURE 3 is an isometric view of the imaging apparatus in FIGURE 2A, 
showing different locations for one or more light source(s); 

FIGURE 4 is an alternative embodiment of the imaging apparatus in 
FIGURES 2A and 3, in which a second set of imaging components and time delay 
5 integration (TDI) detector are included for monitoring light from a bead, to avoid 
interference between different reporters, and showing alternative locations for one or 
more light source(s); 

FIGURE 5 is a schematic diagram illustrating the optical convolution of a 
narrow fluorescence in situ hybridization (FISH) emission spectrum by the present 
1 0 invention, to resolve two FISH probes in a cell; 
Hi FIGURE 6 is a schematic diagram showing the optical convolution of two 

<.n different colors of narrow FISH emission spectra, to resolve the image of the FISH 

H probes on the TDI detector, 

FIGURE 7 is a schematic diagram illustrating how for a wider FISH emission 
Uj 15 spectrum, a deconvolution is provided to resolve the image of two FISH probes of a 
H single color; 

T s FIGURE 8 is a schematic diagram showing the deconvolution of two color 

p FISH spectra that are relatively wide, to resolve the image of the FISH probes; . 

FIGURE 9 is a schematic block diagram of the system used to process the 
; s ;jj 20 signal produced by a TDI detector in the present invention; 

U FIGURE 10 is a schematic diagram illustrating how an imaging system in 

accord with the present invention is used to determine whether a cell is from a male or 
female; 

FIGURE 1 1 is a plan view of an alternate embodiment for an imaging system 
25 usable in the present invention that employs a spectral dispersion component 
comprising a plurality of stacked dichroic filters that spectrally separate the light; 

FIGURE 12 is an X-Y plot of several typical passbands for the dichroic filters 
employed in the embodiment shown in FIGURE 1 1 ; 

FIGURE 13 is a schematic illustration of a detection filter assembly that may 
30 optionally be placed in front of the TDI detector in the embodiment of FIGURE 1 1 to 
. further suppress out-of-band light; 

FIGURES 14A-14E are X-Y plots of transmission vs. wavelength 
corresponding to corresponding passbands of the filter segments of the detection filter 
assembly that may optionally be placed in front of the TDI detector; 
35 FIGURE 15 is a plan view of another embodiment of the configuration of 

FIGURE 11, wherein the spectral dispersion filter system comprises a plurality of 
dichroic cube filters orientated at various angles to create the spectral dispersing effect; 
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FIGURE 16A illustrates another embodiment of the imaging system in accord 
with the present invention, in which the spectral emission is not convolved with the 
image and in which the spectral decomposition occurs in an axis perpendicular to flow 
through the use of separate dichroic filters, imaging lenses, and detectors for each 
5 spectral region; 

FIGURE 16B illustrates yet another embodiment of the imaging system in 
accord with the present invention, in which the spectral emission is not convolved with 
the image and in which the spectral decomposition occurs in an axis perpendicular to 
flow through the use of one imaging lens and separate dichroic Alters, and detectors for 
f9% 10 each spectral region; 

*n FIGURE 16C is an isometric illustration of correction plates added to correct 

Hi for astigmatism induced by a plate beam splitter placed in convergent space; 

FIGURE 17 illustrates images that are projected onto a detector for the spectral 
decomposition embodiment when three beads are in view; 
15 FIGURE 18A is a flow chart illustrating the steps employed in the method of 

the present invention for reading reporter labeled beads in a flow; 

FIGURE 18B is a flow chart illustrating the steps employed in the method of 
the present invention for decoding the reporters on a reporter labeled bead; 

FIGURE 19 is a reporter legend that aids in the identification of reporter 
20 labeled beads; 

FIGURE 20 illustrates hybridizing pligos having sequences that are offset by 
one nucleotide, forming a contiguous overlapping string of sequences; 

FIGURE 21 illustrates the application of the flow imaging sequencing method 
to determine the identity of SNPs; 
25 FIGURE 22 is a flow chart illustrating the steps employed in the method of the 

present invention for constructing sequence contigs for identifying a DNA sequence, 
polymorphic alleles, or expressed genes; 

FIGURE 23 is a plot of channel spectral passbands superimposed over 
emission spectra for common fluorochromes used in reporter beads, carriers, or binding 
30 signals; 

FIGURE 24 is a block diagram illustrating the general steps employed in the 
method of the present invention for applying crosstalk corrections to reporter labeled 
bead imagery prior to decoding the reporter labeled bead assemblies; 

FIGURE 25 is a block diagram for the steps of applying spectral and spatial 
35 corrections to reporter labeled bead imagery; 

FIGURE 26 is a block diagram of the steps of generating the accumulated 
average conelograms between the data images and the reference; and 
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FIGURE 27 is an example of imagery before and after spectral and spatial 
correction. 

Description of the Preferred Embodiment 

It is contemplated that the present invention may be applied to combinatorially 
5 created beads and compounds as well as specifically directed synthesis of beads and 
compounds. Further details regarding both of these aspects of the present invention are 
discussed below. 

New Method for Analyzing Reporter labeled Beads 

A new bead imaging system and method for analyzing reporter labeled beads 

10 addresses the problems for carrying out this task that would arise in reading beads using 
any conventional approach and adds new capabilities to the analysis and handling of 
reporter labeled beads. The imaging system enables the discrimination of the different 
reporters attached to each bead and therefore enables the decoding of the complete 
reporter signature and corresponding chemical identity (e.g., oligonucleotide sequence) 

15 bound to the bead. By preferably handling beads in suspension and using 
hydrodynamic focussing, billion-count bead samples can easily be moved through the 
field of view (FOV) of a flow imaging system at high rates. Existing non-imaging flow 
cytometers cannot be used with this combinatorial scheme because on a given bead, 
different reporters with any colors in common at the same intensity cannot be 

20 distinguished from each other without spatial information. Similarly, without image 
information, different reporter sizes or shapes, which can confer part of the information 
included in a reporter signature, cannot be discriminated. It should be noted that wHile 
a preferred embodiment of the present invention contemplates imaging encoded beads 
while entrained in a flow of fluid, encoded beads can be imaged in stasis, and a flow is 

25 not required. It is contemplated that imaging encoded beads in a flow of fluid will 
enable a large number of beads to be rapidly imaged and analyzed. It is anticipated that 
the resolution afforded by the present invention will enable the reporters themselves to 
be small beads of polystyrene, silica, resin, or any another substance that can be 
fabricated or treated to possess a unique optical signature, which are then associated 

30 with larger (yet still relatively small) substrate beads. 

An embodiment of the present invention for analyzing beads in flow includes 
subsystems for implementing the tasks of: (1) optical signal collection and spectral 
decomposition; (2)pixelated detection; (3) illumination; (4) bead velocity detection; 
and (5) sample handling. One or more embodiments of each of these subsystems is 

35 described in detail below. 

An important aspect of the present invention lies in its ability to simultaneously 
discriminate the location of the various fluorescent emission spectra produced by 
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reporters attached to the solid substrates (beads). This ability enables the rapid 
determination of the presence or absence of a binding event and the identity of each 
bead (i.e., its reporter signature) after a bead has passed through the field of view 
(FOV). Hydrodynamic focussing ensures that the beads are at or near the focal plane of 
5 the imaging system. When required, an optical correction can be applied to account for 
minor changes in the location of the hydrodynamically focussed bead column, hi 
addition, the imaging system can be constructed to view the beads from multiple angles 
to produce image data of a different fraction of the bead surface, to enable the 
construction of a three-dimensional representation of each bead This configuration 
10 also enables correction of focus errors. The present invention spectrally decomposes 
S3 the signal from each reporter in the axis perpendicular to flow and then forms an image 

of the reporters on a bead onto a single detector or multiple detectors. FIGURE 1 
\| schematically illustrates this process. 

As shown in FIGURE 1, a plurality of marked beads 11 pass through a 
15 sampling column 15 in a single file orientation. The imaging system shown includes a 
plurality of detectors 13, focussing elements or lenses 17, and one or more light 
source(s) 19. The location of the image on the detector is then determined by the 
spectra] content of the signal emitted from the reporter as well as the spatial position of 
the reporter with respect to the bead. There are five embodiments of the present 
- - 20 invention that accomplish the spectral decomposition and imaging of the beads. 
First Embodiment for Spectral Decomposition qnd Imaging 

A first embodiment of the system for spectral decomposition and imaging is 
shown in FIGURE 2 A, and this embodiment is also disclosed in commonly assigned 
U.S. Patent No. 6,211,955, entitled 'Imaging and Analyzing Parameters of Small 
25 Moving Objects Such as Cells," filed on March 29, 2000, the drawings and disclosure 
of which are hereby specifically incorporated herein by reference. In this previously 
filed application, there is no discussion of identifying reporters on beads. However, the 
following discussion describes how the apparatus disclosed in this previously filed 
application can be employed for spectral decomposition and imaging of objects such as 
30 beads and reporters included thereon. In regard to the present invention, a bead 
provided with reporters is simply a specific type of object. (It should be noted that 
where there is any variance between the description in any document incorporated 
herein by reference and the present disclosure, the present disclosure takes precedence.) 
Although several types of pixelated detectors may be used in the present 
35 invention, in many cases a time delay integration (TDl) detector is preferable. The TDI 
detector preferably comprises a rectangular charge-coupled device (CCD) that employs 
a specialized pixel readout algorithm, as explained below. Non-TDI CCD arrays are 
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commonly used for two-dimensional imaging in cameras* In a standard CCD array, 
photons that are incident on a pixel produce charges thai are trapped in the pixel The 
photon charges from each pixel are readout of the detector array by shifting the charges 
from one pixel to the next, and then onto an output capacitor, producing a voltage 
5 proportional to the charge. Between pixel readings, the capacitor is discharged and the 
process is repeated for every pixel on the chip. During the readout, die array must be 
shielded from any light exposure to prevent charge generation in the pixels that have 
not yet been read. 

In one type of TDI detector 44, which comprises a CCD array, the CCD array 
10 remains exposed to the light as the pixels are readout The readout occurs one row at a 
time from the top toward the bottom of the array. Once a first row is readout, the 
remaining rows are shifted by one pixel in the direction of toe row that has just been 
read. If the object being imaged onto the array moves in synchrony with the motion of 
fji the pixels, light from the object is integrated for the duration of the TDI detector's total 

1 5 readout period without image blurring. The signal strength produced by a TDI detector 
will increase linearly with the integration period, which is proportional to the number of 
TDI rows, but the noise will increase only as the square root of the integration period, 
resulting in an overall increase in the signal-to-noise ratio by the square root of the 
number of rows. One TDI detector suitable for use in the present invention is a Dalsa 
20 Corp., Type IL-E2 image sensor, although other equivalent or better image sensors can 
alternatively be used. The Dalsa image sensor has 96 stages or rows, each comprising 
512 pixels; other types of image sensors useable in the present invention may have 
different configurations of rows and columns or a non-rectilinear arrangement of pixels. 
The Dalsa sensor has approximately 96 times the sensitivity and nearly 10 times the 
25 signal-to-noise ratio of a standard CCD array. The extended integration time associated 
with TDI detection also serves to average out temporal and spatial illumination 
variations, increasing measurement consistency. 

It should be emphasized that the present invention is not limited to TDI 
detectors or CCD types of TDI detectors. Other types of TDI detectors, such as 
30 complementary metal oxide semiconductor (CMOS) and multi-channel plate imaging 
devices might also be used for the TDI detector in the present invention. It is important 
to understand that any pixela&d device (i.e., having a multitude of light sensitive 
regions) in which a signal produced in response to radiation directed at the device can 
be caused to move through the device in a controlled fashion is suitable for use as the 
35 TDI detector in the present invention. Typically, the signal will move in synchrony 
with a moving image projected onto the device, thereby increasing the integration time 
for the image, without causing blurring. However, the motion of the signal can be 
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selectively desynchronized from the motion of the radiation image, as required to 
achieve a desired affect. 

In imaging systems 45, 70, 250 in FIGURES 2C, 4, and 1 1 respectively, and in 
other embodiments of the present invention that employ a fluid flow to carry objects 
through the imaging system, a flow-through cuvette or a jet (not shown) contains the 
reporter labeled beads or other objects being analyzed The velocity and cellular 
concentration of the fluid may be controlled using syringe pumps, gas pressure, or other 
pumping methods (not shown) to drive a sample solution through the system to match 
the pixel readout rate of the TDI detector. However, it should be understood that the 
readout rate of the TDI detector can be selectively controlled, as required, to match the 
motion of the sample solution. 

As will be evident in FIGURE 2B, if the Figure depicts the imaging of 
object 24 over time, the object is shown at both a position 26 and a position 28 as it 
moves with fluid flow 22. As a consequence, images of object 24 will be produced on 
the detector at two discrete spatial positions 26' and 28', as indicated on the right side of 
FIGURE 2B. In the case of TDI detection, the signal on the detector can be made to 
move in synchronicity with image of object 24 such that signal is collected and 
integrated over the entire traversal of the detector. Alternatively, if FIGURE 2B is 
depicting a single instant in time, positions 26 and 28 can represent the location of two 
separate objects, which are simultaneously imaged on the detector at positions 26' 
and 28'. 

As shown in the FIGURE 2A, light from a column 22 of beads 24 that is 
hydrodynamicaUy focussed to a well defined region, is collimated by passing through a 
collection lens 32. The light from the beads travels along a collection path 30. A 
spectral dispersing element 36 disposed in the collection path spectrally disperses the 
collimated light that has passed through the collection lens in a plane substantially 
orthogonal to a direction of relative movement between the beads and the imaging 
system, producing spectrally dispersed light An imaging lens 40 is disposed to receive 
the spectrally dispersed light, producing an image from the spectrally dispersed 
light 38. Also included is a pixelated detector 44, disposed to receive the image 
produced by the imaging lens. As the movement occurs, the image of the bead 
produced by the imaging lens moves from row to row across the detector! As will be 
described earlier, the detector may be a TDI-type detector or a frame type detector, 
depending upon the specific embodiment. 

As a result of light collimation by the collection lens in this embodiment, all 
light emitted from a first point in the bead travels in parallel rays. U ght emitted from a 
second point in the bead will also travel in parallel rays, but at a different angle relative 
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to light from the first point. In this manner, spatial information in the bead is 
transformed by the collection lens into angular information in the collection path. The 
spectral dispersing element acts on the collimated light such that different spectral 
components leave the spectral dispersing element at different angles, in a plane 
5 substantially orthogonal to the direction of the relative movement between the bead and 
the imaging system. In this manner, both spatial and spectral information in the bead 
are transformed into angular information. The imaging lens acts on the light from the 
dispersing element to transform different light angles into different positions on the 
detector. Spatial information is preserved by the system since light from the different 
10 positions in the bead is projected along a path 42 to different positions on the detector, 
in both axes. In addition, light of different spectral composition that originates from the 
bead is projected to different positions on the detector in an axis substantially 
!* orthogonal to the relative movement between the bead and the imaging system. In this 

fj\ manner, the spatial information from the bead is preserved, and spectral information 

f « 1 5 covering a large bandwidth at high resolution is simultaneously collected. 
v{ When used for bead identification in accord with the present invention, this 

[t apparatus provides substantial utility in resolving reporter locations and spectra on the 

^ detector, even when the reporters are disposed in spatially close relationship within a 

bead. When spectral imaging occurs in the present invention, the spatial distribution of 
rij 20 light in the bead is convolved with the spectral distribution of that light to produce the 
^5 image of the bead at the detector. This convolution can result in blurring in the 

^ dispersion axis, depending on the spectral bandwidth of the light Narrow spectral 

bandwidths will result in little or no blurring, depending on the spectral resolution of the 
system. Li the present invention, it is contemplated that the spectral resolution will be 
25 approximately 3 nm per pixel, with a spatial resolution in object space of approximately 
1 micron. However, the spatial and spectral resolution can be adjusted to match the 
requirements of a particular application and the values noted herein should not be 
considered as limiting on the scope of the present invention. 

Alternate embodiments of the imaging system illustrated in FIGURE 2A and 
30 useful for imaging beads are shown in FIGURES 2G4. These alternate embodiments 
differ in the number and orientation of various optical components, as described below, 
but generally function in the same manner as the imaging system of FIGURE 2A. 

FIGURE 2C illustrates an imaging system 45 that includes a slit 52 used to 
prevent extraneous light from reaching pixelaled detector 44. In imaging system 45, 
35 light 46 from bead 24 (or other object) is focussed by an objective lens 48 onto slit 52. 
Slit 52, as shown in FIGURE 2C, is sufficiently narrow to block light that is not 
focussed onto the slit by objective lens 48, thereby preventing extraneous light from 
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passing through the slit Ligftt30', which has passed through the slit, is collected by 
collection lens 32, as discussed above, producing collected light 34. Collected light 34 
is spectrally dispersed by prism 36 and is imaged by imaging lens 40 onto pixelated 
detector 44, also as discussed above. By excluding light other than that from object 24 
from reaching pixelated detector 44, the detector produces an output signal that 
corresponds only to the actual images of an encoded bead, and the signal is not affected 
by the extraneous light that has been excluded. If not excluded in this manner, ambient 
ligfrt reaching pixelated detector 44 might otherwise produce **noise" in the output 
signal from the pixelated detector. 

It should be noted that in the illustration of each of the imaging systems in 
FIGURES 2A and 2C, a light source has not been shown. These firet two 
embodiments have been illustrated in their most general form to make clear that a 
separate light source is not required to produce an image of an encoded bead if the 
object is luminescent, i.e„ if the object produces light. However, in most uses of the 
imaging apparatus for imaging beads, one , or more light source(s) will be used to 
provide light that is incident on the bead being imaged. The location of the light 
sources substantially affects the interaction of the incident light with the bead (or other 
object) and the kind of infomiation that can be obtained from the images on the 
pixelated detector. 

In FIGURE 3, several different locations of light sources usable to provide light 
incident on bead 24 are illustrated It should be understood, however, that light sources 
can be located at many other positions besides those shown in FIGURE 3. The location 
of each light source that is employed will be dependent upon the kind of imaging of the 
bead, and the kind of data for the bead, to be derived from the signal produced by the 
pixelated detector. For example, employing a light source 60a or a light source 60b, as 
shown in the figure, will provide light 58 that is incident on bead 24 and which is 
scattered from the bead into the optical axis of collection lens 32. The optical axis of 
collection lens 32 is at about a 90° angle relative to the directions of the light incident 
upon bead 24 from either light source 60a or 60b. In contrast, a light source 62 is 
disposed so that light 58 emitted from that source travels toward the bead in a direction 
that is generally aligned with the optical axis of collection lens 32, so that the image 
formed on detector 44 will not include light absorbed by bead 24. Light absorption 
characteristics of the bead can thus be determined by illuminating the bead using light 
source 62. 

A light source 64 is disposed to illuminate bead 24 with light directed toward 
the bead along a path that is approximately 3045° off the optical axis of collection 
lens 32. This light 58, when incident on bead 24, will be reflected (scattered) from 
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bead24, and the reflected or scattered light will be imaged on detector 44. A more 
directly reflected light is provided by an epi light source 66, disposed so as to direct its 
light 58 toward a partially reflective surface 68 that is disposed so that a portion of the 
light is reflected through collection lens 32 and onto bead 24. The light reaching the 
bead will be reflected from it back along the axis of collection lens 32 and will at least 
in part pass through partially reflective surface 68 to form an image of the bead on 
pixelated detector 44. Alternatively, a dichroic mirror may be employed instead of, and 
in the position of, partially reflective surface 68 to direct light from epi light source 66 
to excite fluorescence or other stimulated emission ftom bead 24. Light emitted from 
bead 24 at a different wavelength than the light source is then at least partially collected 
by collection lens 32 and passes through the dichroic minor for spectral dispersion and 
detection by the detector. 

Each of the light sources illustrated in FIGURE 3 produce light 58, which can 
either be coherent, noncoherent, broadband, or narrowband light, depending upon the 
application of die imaging system for imaging beads. Thus, a tungsten filament light 
source can be used for applications in which a narrowband light source is not required. 
For applications such as stimulating the emission of fluorescence ftom reporters, 
narrowband light is preferred, since it also enables a spectrally-decomposed, 
nondistorted image of the bead to be produced from light scattered by the bead This 
scattered light image will be separately resolved from the reporters produced on 
pixelated detector 44, so long as the emission spectra of any reporters are at different 
wavelengths than the wavelength of the light The light source can be either of the 
continuous wave (CW) or pulsed type. If a pulsed type illumination source is 
employed, the extended integration period associated with pixelated detection can allow 
the integration of signal from multiple pulses. Furthermore, it is not necessary for the 
light to be pulsed in synchronization with the pixelated detector. 

FIGURE 4 illustrates an arrangement that enables the imaging of a bead from 
two different directions, in order to distinguish features that would otherwise overlap 
when viewed from a single direction. A stereoscopic imaging system 70 in FIGURE 4 
includes two pixelated detectors 44a and 44b, and their associated optical components, 
as discussed above in connection with the imaging system of FIGURE 2A. 

By positioning the optical axes of collection lenses 32 for the two pixelated 
detectors so that they are spaced apart, for example, by 90°, it is possible to separately 
resolve the signature of reporter tags imaged from two or more reporters on at least one 
of pixelated detectors 44a or 44b. If two or more reporters overlap in regard to the 
image produced on one of the detectors, they will be separately resolved in spectrally 
dispersed image produced on the other pixelated detector. Further, the use of two 
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pixelated detectors in imaging system 70 in what might be referred to as a "stereo or 
three-dimensional configuration" allows flexibility in the configuration of each leg of 
the system, including parameters such as the relative pixelated readout rates, axial 
orientations, inclinations, focal plane positions, and magnification. Multiple beads or 
other objects may be imaged onto each detector simultaneously in the vertical direction. 
Since the beads may move in synchronicity with the signal on the pixelated detector, no 
gate or shutter is required to prevent blurring of the image. As previously noted, the 
present invention can use a pulsed or CW light source without need for a trigger 
mechanism to time a pulse coincident with particle arrival in the FOV, 

Also illustrated in FIGURE 4 are several exemplary positions for light sources, 
which are useful for different purposes in connection with the imaging system 
illustrated therein. In connection with pixelated detector 44a, light source 62 provides 
illumination of bead 24 from a direction so that absorption characteristics of the bead 
can be determined from the image produced on the pixelated detector. At the same 
time, light provided by light source 62 that is scattered from bead 24 can be used to 
produce a scatter image and spectrally dispersed images on pixelated detector 44b. A 
light source 74 can be employed to produce spectrally dispersed and scattered images 
on both pixelated detectors 44a and 44b. If light sources 62 and 72 are of different 
wavelengths and an appropriate filter is provided to block the wavelength from the light 
source aligned with the optical axis of the respective collections lenses 32, these two 
light sources can be used for producing scattered light from the bead. For example, 
suppose light source 72 produces light of a wavelength A that scatters ftom bead 24 
and is directed toward pixelated detector 44a By including a filter (not shown) that 
blocks wavelength B produced by light source 62, the light at wavelength B will not 
directly affect the images produced on pixelated detector 44a. Similarly, the light from 
light source 72 would be blocked with an appropriate filter (not shown), so that it does 
not interfere with the imaging of light produced by light source 62 that is scattered from 
bead 24 onto pixelated detector 44b. 

Epi light source 66 is also illustrated for use in producing images on pixelated 
detector 44a in conjunction with partial reflector 68. light source 64 can be used to 
generate reflected light to produce images on pixelated detector 44a, while scattered 
light from this source is directed toward pixelated detector 44b, These and other 
possible locations of light sources will be apparent to those of ordinary skill in the art, 
as appropriate, for providing the incident light on the bead needed to achieve imaging, 
depending upon the particular application and information about the bead that is 
desired. 
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FIGURE 5 illustrates a detector used in the present invention, with a spectral 
resolution of approximately 10 nm per pixel and a spatial resolution of approximately 
0.5 microns, hi the following discussion of FIGURES 5-8, the operation of the 
imaging system is directed toward the identification of fluorescence in situ 
hybridization (FISH) probes bound to specific DNA with cells. However, those skilled 
in the art will appreciate that the same method and apparatus applies to the spectral and 
spatial information resulting from the imaging of reporters associated with beads. With 
respect to the following discussion, a cell and nucleus can be considered to be 
equivalent to a bead, and FISH probes can be considered to be equivalent to reporters 
associated with or bound to a bead. 

FIGURE 5 illustrates how die imaging system that is usable to image reporter 
beads is used to image a cell 140 having a nucleus 142 in which are disposed two FISH 
probes 144a and 144b having the same emission spectrum. As noted above, in a bead 
analysis, the FISH spots may be considered analogous to reporters, m FIGURE 5, the 
emission spectrum 146 of FISH probes 144a and 144b is approximately 10 nm in 
width, such as would be produced by "quantum dots" or a narrowband fluorescent dye. 
The optical convolution of the narrow bandwidth spectrum results in minimal blurring 
of FISH spots 148a and 148b, enabling them to be readily resolved on detector 44. 

m FIGURE 6, a cell 150 is iOustrated having a nucleus 152 in which are 
disposed FISH probes 154 and 156 having different emission spectra. Each of the 
emission spectra of FISH probes 154 and 156 are relatively narrow, such as the 
emission spectra mom quantum dots, as indicated by wavebands 158 and 160, and 
therefore, as in FIGURE 5, minimal blurring occurs in FISH spots 162 and 164. 
Furthermore, the spectral dispersion of the present invention, which maps wavelength 
into lateral position on TDI detector 44, produces a relatively wide physical 
displacement of FISH spots 162 and 164, despite the close proximity of FISH 
probes 154 and 156 in the cell. Taken together, FIGURES 5 and 6 illustrate how the 
present invention discriminates FISH probes of the same or different color, thereby 
enabling the simultaneous enumeration of numerous genetic traits. FIGURE 10 
illustrates how a cell and FISH spots contained therein will be imaged by the present 
invention upon a pixelated detector. 

FIGURES 7 and 8 illustrate that the present invention can also be used with 
light of wider spectral bandwidth. In this case, an additional signal processing step is 
performed to correct for lateral blurring due to the wide emission spectra. Note that this 
involves a deconvolution step which requires the use of a pixelated detector that is a 
TDI detector. In FIGURE 7, a cell 140 having a nucleus 142 is shown, and FISH 
probes 170a and 170b having a common emission spectrum are disposed in the 
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nucleus. FISH probes 170a and 170b are characterized by producing a relatively wide 
emission spectrum 172. When optically convolved by the spectral dispersion provided 
by the present invention, FISH spots 174a and 174b are produced on a TDI 
detector 44a, but their images are laterally blurred across pixelated detector 44, as a 
result of their relatively wide emission spectrum. To more clearly resolve the 
separation of FISH spots probes 174a and 174b, a deconvolution is carried out on the 
signal produced by TDI detector 44a, with the known FISH emission spectrum, thereby ' 
producing accurate FISH spot representations 178a and 178b on a display 176. The 
deconvolution step enhances the ability to enumerate the number of FISH spots. Note 
that the deconvolution step requires the use of a pixelated detector that is a TDI 
detector. 

FIGURE 8 illustrates a corresponding relationship between FISH probes 180 
and 182, which are disposed within a nucleus 152 of a cell 150. FISH probes 180 
and 182 are characterized by each producing relatively wide band emission spectra 184 
and 186, as shown in the Figure. Optical convolution of the fluorescence emitted by 
the FISH probes, which are spectrally dispersed, produces FISH spots 188 and 190 on 
TDI detector 44. Again, by deconvolving the known FISH emission spectra with the 
signal produced by TDI detector 44, the corresponding images shown on display 176 of 
FISH spots 192 and 194 are recovered. Again, the spectral dispersion of the present 
invention, which maps wavelength into lateral position on TDI detector 44, produces a 
relatively wide physical displacement of FISH spots 192 and 194, despite the close 
proximity of FISH probes 180 and 182 in the cell. In this manner, it is possible to 
resolve these images of FISH spots produced by FISH probes having different and 
relatively wide emission spectra. 

A system 230 for analyzing the signal produced by TDI detector 44a and 
performing the deconvolution steps described above is illustrated in FIGURE 9. In this 
figure, the signal from TDI detector 44a is applied to an amplifier 232, which buffers 
the signal and amplifies it to achieve a level required by an analog-to-digital (A/D) 
converter 234. This A/D converter converts the analog signal from amplifier 232 into a 
digital signal that is input into a line buffer 236. line buffer 236 temporarily stores the 
digital signal until it can be processed by a CPU 238. To cany out the deconvolution 
noted above, a spectral buffer 240 is loaded with the known emission spectrum for each 
of the FISH probes being used so that their emission spectra can be deconvolved with 
the signal stored in line buffer 236. CPU 238 is a high speed processor programmed to 
carry out the deconvolution and other analysis procedures, enabling the identification of 
desired characteristics or parameters of the object being imaged. The output from 


flIOIi»l7.|-l2dfOOI7AP.<l« 


-22- 


CPU 238 is temporarily stored in an image line buffer 242 that enables the image to be 
displayed or otherwise recorded for later analysis. 

FIGURE 10 illustrates a practical application of the above described systems 
for identifying a male cell 200 and a female cell 208 and for producing their 
5 corresponding scatter images 212 and 220. Male cell 200 includes a nucleus 202 that 
has been stained with a yellow fluorescent dye. In addition, a FISH probe 204 
produces a fluorescent orange emission, indicating the presence of an X-chromosome 
in the nucleus, while a FISH probe 206 produces red fluorescence emission, indicating 
the presence of a Y-chromosome. Spectral decomposition of the fluorescence 
10 emissions from male cell 200, when the cell is illuminated with light from a green light 
source, results in a series of images on TDI detector 44, separated as a function of the 
wavelength of the light that is imaged Green light that is incident on the cells has a 
narrow waveband, and image 212 of male cell 200 produced by green light scatter is 
only slightly convoluted by the spectral decomposition process. Green light scatter 
UJ 15 image 212 of cell 200 and its nucleus 202 appear on the left side of the TDI detector, 
J :| while a fluorescent spot 214 corresponding to the yellow fluorescence emitted by 

ji* nucleus 202 appears in the next few columns on the TDI detector. Furthermore, as a 

0 function of the different wavelengths of the fluorescence emitted by FESH probes 204 

|;y and 206, FISH spots 216 and 218 appear at locations spaced apart on the detector, but 

q 20 slightly blurred across the columns of TDI detector 44 due to the widths of their 
respective emission spectra. By analyzing the signals produced by the TDI detector, 
the FISH probes responsive to X and Y chromosomes are detected, enabling the user to 
determine that cell 200 is a male cell, since it includes both the X and Y chromosome. 
Similarly, female cell 208, when spectrally decomposed, also includes the characteristic 
25 yellow fluorescence of nucleus 2 10, but unlike the male cell, includes two FISH 
spots 216 corresponding to FISH probes 204, which indicates the presence of two X- 
chromosomes. Because TDI detector 44 also distinguishes the spatial position of male 
cell 200 and female cell 208, the corresponding spectral decompositions for these cells 
are readily separately resolved as both cells pass through the imaging system in the 
30 direction indicated by the arrow to the lower left of FIGURE 10. Again, it should be 
noted that a deconvolution can be applied to the signal produced by TDI detector 44 to 
provide better resolution of the corresponding FISH spots that are illustrated 

Again, to avoid any confusion, those skilled in the art will appreciate that in 
the present invention, which is employed for the analysis of reporter labeled beads, a 
35 bead in analogous to the cell and a reporter is analogous to a FISH spot, in regard to 
the discussion of FISH spots contained herein. One difference is that a FISH spot is 
typically constructed with a single fluorochrome and therefore emits over a single 
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spectral range, as will be evident in the figures. In the case of reporters, each reporter 
may be constructed of more than one fluorochrome and therefore, can emit several 
spectra simultaneously. The same signal processing techniques disclosed in the 
above-referenced U.S. Patent No. 6,21 1,955 and discussed relative to objects and 
5 FISH spots can be applied to enumerate unique reporters of more than one 
fluorochrome associated with, or bound to a bead. Therefore, it will be understood 
that the imaging system and techniques used for analysis may be used to enumerate 
unique reporters associated with or bound to a bead in order identify the unique 
spectral signature of the bead 
10 FIGURE 17 illustrates the results of processing the images on the detector in a 

reporter labeled bead scenario, wherein the emission spectra from each reporter has 
been deconvolved. FIGURE 17, which is discussed in greater detail below, is also used 
to illustrate the use of the non-convolving spectral decomposition and imaging system 
discussed in the next section. 
Ul 15 Second and Third Embodiments of Apparatus for Spectral Decomposition and Imaging 

t, 7 

In the second embodiment of apparatus usable in practicing the present 
invention, a spectral dispersing component having characteristics that ensure no 
distortion or convolution of the image occurs due to the emission bandwidth is 
employed, and as a result, a deconvolution step is not needed to process the image data, 
i 20 FIGURE 1 1 illustrates this second preferred embodiment Unlike a prism, where every 
wavelength leaves the prism at a different angle, all light within a predefined bandwidth 
incident on the dichroic beam splitter at a common angle leaves the dichroic beam 
splitter at the same angle. Consequently, there is no convolution between the emission 
spectrum of the light leaving the bead and the image of that bead. When using such a 
25 spectral dispersing component, light of a first spectral bandwidth is reflected from the 
first dichroic beam splitter at a predefined nominal angle. light of a second spectral 
bandwidth is passed through the first dichroic beam splitter to the next dichroic beam 
splitter and is reflected therefrom at a different predefined nominal angle, light of a 
third spectral bandwidth is passed through the first and second dichroic beam splitters 
30 to a third dichroic beam splitter and reflected therefrom at a third predefined nominal 
angle. The dichroic beam splitters are selected to cover the desired light spectrum with 
the appropriate spectral passbands. In FIGURE 11 illustrates a five-color stacked 
wedge spectral dispersing filter assembly 252. This second embodiment is 
substantially similar to the embodiment shown in FIGURES 2 and 3, except that 
35 spectral dispersing prism element 36 (of FIGURES 2 and 3) is replaced by spectral 
dispersing filter assembly 252. The spectral dispersing filter assembly splits the light 
into a plurality of light beams having different band widths. Each light beam thus 
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produced is directed at a different nominal angle so as to fall upon a different region of 
detector 44. The nominal angular separation between each bandwidth produced by the 
spectral dispersing filter assembly 252 exceeds the field angle of the imaging system in 
object space thereby preventing overlap of the field images of various bandwidths on 
the detector. 

Spectral dispersing filter assembly 252 comprises a plurality of stacked dichroic 
wedge filters, including a red dichroic filter R, an orange dichroic filter O, a yellow 
dichroic filter Y, a green dichroic filter G, and a blue dichroic filter B. Red dichroic 
filter R is placed in the path of collected light 34. oriented at an angle of 
approximately 44.0° relative to an optic axis 253 of collection lenses 32 and 32. Light 
of red wavelengths and above, i.e., >640nm, is reflected from the surface of red 
dichroic filter R at a nominal angle of 1°, measured counter-clockwise from a vertical 
optic axis 257. Example spectral reflectance characteristics R' of red dichroic filter R 
are plotted in FIGURE 12, along with example spectral reflectance characteristics 
corresponding to the other dichroic filters used in spectral dispersing filter 
assembly 252. In FIGURE 12, O' indicates the spectral reflectance characteristics of 
orange dichroic filter O, Y indicates the spectral reflectance characteristics of yellow 
dichroic filter Y, etc. The light reflected by red dichroic filter R leaves spectral 
dispersing filter assembly 252 and passes through imaging lenses 40a and 40b, which 
cause the light to be imaged onto a red light receiving region of TDI detector 44, which 
is disposed toward the right end of the TDI detector, as shown in FIGURE 1 1 . 

Orange dichroic filter O is disposed a short distance behind red dichroic filter R 
and is oriented at an angle of 44.5 degrees with respect to optic axis 253. light of 
orange wavelengths and greater, i.e., > 610 nm, is reflected by orange dichroic filter O 
at a nominal angle of 0.5° with respect to vertical optic axis 257. Because the portion 
of collected light 34 comprising wavelengths longer than 640 nm was already reflected 
by red dichroic filter R, the light reflected from the surface of orange dichroic filter O is 
effectively bandpassed in the orange colored region between 610 nm and 640 nm. This 
light travels at a nominal angle of 0.5° from vertical optic axis 257, and is imaged by 
imaging lenses 40a and 40b so as to fall onto an orange light receiving region disposed 
toward the right-hand side of TDI detector 44 between a center region of the TDI 
detector and the red light receiving region. 

Yellow dichroic filter Y is disposed a short distance behind orange dichroic 
filter O and is oriented at an angle of 45° with respect to optic axis 253. Light of yellow 
wavelengths, i.e., 560 nm and longer, is reflected from yellow dichroic filter Y at a 
nominal angle of 0.0° with respect to vertical optic axis 257. Wavelengths of light 
reflected by yellow dichroic filter Y are effectively bandpassed in the yellow region 
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between 560 nm and 610 nm and are imaged by imaging lenses 40a and 40b near 
vertical optic axis 257 so as to fall on a yellow light receiving region toward the center 
of TDI detector 44. 

In a manner similar to dichroic filters R, O, and Y, dichroic filters G and B are 
5 configured and oriented so as to image green and blue light wavebands onto respective 
green and blue light receiving regions of TDI detector 44, which are disposed toward the 
left-hand side of the TDI detector. By stacking the dichroic filters at different predefined 
angles, spectral dispersing filter assembly 252 collectively works to focus light within 
predefined wavebands of the light spectrum onto predefined regions of TDI detector 44. 
10 Those of ordinary skill in the art will appreciate that the filters used in the spectral 
dispersing filter assembly 252 may have spectral characteristics that differ from those 
described above and in FIGURE 12. Further, the spectral characteristics may be 
arbitrary and not limited to dichroic in order to achieve the desired dispersion 
t .» characteristics. 

ry 1 5 The wedgp shape of the dichroic filters in die preceding discussion allows the 

filters to be placed in near contact, in contact, or possibly cemented together to form the 
spectral dispersing filter assembly 252. The angle of the wedge shape fabricated into 
the substrate for the dichroic filter allows easy assembly of the spectral dispersing filter 
assembly 252, forming a monolithic structure in which the wedge-shaped substrate is 
r i j 20 sandwiched between adjacent dichroic filters. If the filters are in contact with each 
Q other or cemented together, the composition of the materials that determine the spectral 

performance of the filter may be different from those which are not in contact. Those 
of ordinary skill in the an will appreciate that flat, non wedge-shaped substrates could 
be used to fabricate the spectral dispersing filter assembly 252. In this case another 
25 means such as mechanically mounting the filters could be used to maintain the angular 
relationships between the filters. 

In addition to the foregoing configuration, non-distorting spectral dispersion 
system 250 may optionally include a detector filter assembly 254 to further attenuate 
undesired signals in each of the light beams, depending upon the amount of rejection 
30 required for out-of-band signals, FIGURE 13 illustrates the construction of an 
exemplary detector filter 254 corresponding to the five color bands discussed above and 
includes a blue spectral region 256, a green spectral region 258, a yellow spectral 
region 260, an orange spectral region 262, and a red spectral region 264, all of which 
are disposed side-by-side, as shown in the Figure. Hie corresponding spectral 
35 characteristics of the blue, green, yellow, orange, and red spectral regions or wavebands 
are respectively shown in FIGURES 14A-14E The detection filter assembly shown in 
FIGURE 13 may be constructed by cementing separate filters in side-by-side 
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arrangement on a common substrate or by other means well known to those or ordinary 
skill in the art Additionally, the ordinary practitioner in the art will understand that the 
filter may alternatively be placed at an intermediate image plane, instead of direcdy in 
front of TDI detector 44. 

In the embodiment shown in FIGURE 1 1, light may pass through each dichroic 
fflter in the spectral dispersing filter assembly 252 twice before exiting the spectral 
dispersing filter assembly 252. This condition will further attenuate out-of-band 
signals, but will also attenuate in-band signals. FIGURE 15 illustrates a third 
embodiment useful in the present invention in which the light does not pass through 
another dichroic filter after reflection, hi this embodiment, a plurality of cube dichroic 
filters, including a red cube filter 266, a yellow cube filter 268, a green cube filter 270, 
and a blue cube filter 272 are spaced apart sufficiendy to ensure that light does not pass 
through any of the cube filters more than once. As with the second embodiment of 
FIGURE 1 1, the cube dichroic filters are oriented at appropriate angles to image light 
within a predefined bandwidth to distinct regions on a TDI detector 274. As the light is 
reflected from each of cube dichroic filters 266, 268, 270 and 272, it is directed toward 
imaging lenses 40a and 40b, and different bandpass portions of the light are focussed 
upon corresponding red, yellow, green, and blue light* receiving segments or regions 
defined on a light receiving surface of TDI detector 274. If desired, an optional 
detector filter assembly 276 of similar construction to detector filter assembly 254 (but 
without the orange spectral region) may be used to increase the rejection of out-of-band 
signals, It should be apparent to those skilled in the art that separate spaced apart plates, 
or pellicle beam splitters could also be used in this application instead of the cube 
filters. 

In the third embodiment illustrated in FIGURE 15, the image lenses 40a 
and 40b must be placed a sufficient distance away from the plurality of cube filters to 
minimize the clear aperture requirement for lenses 40a and 40b. Those skilled in the 
art will appreciate the clear aperture in the plane orthogonal to the page can increase as 
the distance between the lenses and plurality cube filters increases. Therefore, the 
placement of lenses 40a and 40b must be chosen to appropriately accommodate the 
clear aperture in both planes. 

The foregoing descriptions of the preceding non-convolving embodiments 
illustrate the use of four and five color systems. Those skilled in the art will appreciate 
that a spectral dispersing component with more or fewer filters may be used in these 
configurations in order to construct a system covering a wider or a narrower spectral 
region, or different passbands within a given spectral region. Likewise, those skilled in 
the art will appreciate that the spectral resolution of the present invention may be 
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increased or decreased by appropriately choosing the number and spectral 
characteristics of the dichroic and or bandpass filters that are used. Furthermore, those 
skilled in the art will appreciate that the angles or orientation of the filters may be 
adjusted to direct light of a given bandwidth onto any desired point on the TDI detector 
5 In addition, there is no need to focus the light in increasing or decreasing oider by 
wavelength. For example, in fluorescence imaging applications, one may wish to 
create more spatial separation on the TDI detector between the excitation and emission 
wavelengths by changing the angles at which the filters corresponding to those 
wavelengths are oriented with respect to the optic axes of the system. Finally, it will be 
* n 10 clear to those skilled in the art that dispersion of the collected light may be performed 
on the basis of nonspectral characteristics, including angle, position, polarization, phase, 
or other optical properties. 
{ ^ FIGURE 17 illustrates the images projected onto a detector for the present 

f y spectral decomposition embodiment in the case where three beads are in view. In this 

15 illustration, each bead has a series of four unique reporters visible in which the reporters 
, " are constructed of up to four distinct fluorochromes. There is no spreading of the 

h= reporter images, as occurs in the prism based spectral decomposition embodiment. The 

field angle orthogonal to flow in object space is also indicated on FIGURE 17. In this 
1*11 particular configuration, the field angle in object space is less than +/- 0.25°. Those 

O 20 skilled in the art will appreciate that the field angle can be made larger or smaller. To 
^ the extent that the field angle is made larger, for example, to image cells over a wider 

region on a slide or in a broad flat flow, the field angle at the detector will increase in 
proportion to the number of colors used. Those skilled in the art will appreciate that 
broad flat flow can easily be created using commercially available flow cells as shown 
25 in FIGURE 16B containing flow cell 290, with a cross section that is elongated in an 
axis perpendicular to both the flow and optical axes. The generation of broad flat flow 
is discussed in many references including U.S. Patent No, 5,422,712. 
Fourth Embodiment of Apparatus for Spectral Decomposition and Imaging 

The fourth embodiment of apparatus usable for spectral decomposition and 
30 imaging in connection with the present invention, which is illustrated in FIGURE 16A, 
is similar to the second and third embodiments in that no convolution of the emission 
spectra with the image occurs as a result of the spectral decomposition process. 
Spectral decomposition occurs in an axis generally perpendicular to flow through the 
use of dichroic filters as previously described. However, in this embodiment, separate 
35 imaging lenses and detectors are used for each spectral region. Dichroic filters 277 are 
placed at infinity (with respect to the object) after a collection lens 278 to minimize 
optical aberrations. After each dichroic filter 276a separate imaging lens 280 is used to 
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form an image of the objects onto the corresponding detectors 282. In this 
configuration, each detector 282 has fewer pixels than otter embodiments because each 
detector covers only one spectral band. In the case where a six-color version of this 
embodiment is employed, the images on a single detector would appear like the images 
5 seen on one zone of the detector illustrated in FIGURE 17. For example, the images 
seen on the detector configured to receive light in the red part of the spectrum would 
appear like the right-most zone of FIGURE 17. Since the total number of pixels in 
each detector is low, these detectors may operate at very high speeds. The embodiment 
shown in FIGURE 16A has the advantage of being more efficient than other 

10 embodiments, because the light from the object only passes through each dichroic 
filter 277 once. A still further advantage of this embodiment is that each detector 282 
can be focused independently for each color thereby simplifying the optical design by 
removing the constraint of longitudinal color correction. An additional objective 
lens 48 and slit 52, such as that shown in FIGURE 2C can be incorporated into the 

15 embodiment shown in FIGURE 16A to prevent unwanted light from striking one or all 
of the pixelated detectors. 

Fifth Embodiment of Apparatus for Spectral Decomposition and Imaging 

FIGURE 16B illustrates another embodiment for spectral decomposition 
imaging. It is similar to the previous embodiment, however, it has the advantage of 

20 reducing the number of imaging lenses required to project an image upon the detectors. 
In the imaging system of FIGURE 16B, light from beads (or other objects) is focused 
by an objective lens 281 onto a slit 283. Slit 283 is sufficiently narrow to block light 
that is not focussed onto the slit by objective lens 281, thereby preventing extraneous 
light from passing through the slit light which has passed through the slit, is collected 

25 by a collection lens 284 and then directed to imaging lens 285, berth of which are placed 
before the dichroic filters 286. Each detector 287 is placed at the appropriate position 
along the optical path to image the object onto the surface of each detector 287. 
Detectors 287 are placed at the appropriate positions along the optical path to image 
beads or other objects onto detectors 287. The filters are placed in convergent space 

30 with respect to the image of the object and therefore each filter, depending upon its 
design, may impart astigmatism, coma, spherical and chromatic aberration into the 
imagery at each downstream detector. Progressively more of each of these aberrations 
are added by each subsequent filter. In a typical implementation of the present 
invention, the numerical aperture in the filter space is approximately 0.03. Therefore, if 

35 cube substrates are employed for the dichroic filters, coma and astigmatism are 
negligible while spherical aberration is less than 0.15 waves peak. Longitudinal 
chromatic aberration is effectively canceled by moving the detectors to the plane of best 
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focus for their respective color band. Pellicles can also be used in place of the cubes 
with excellent theoretical optical performance. 

If plate substrates are employed for dichroic filters 286, then astigmatism 
dominates the aberrations. Astigmatism is imparted in the transmitted wavefront 
5 through a plate dichroic filter, but can effectively be cancelled by inserting a clear 
correction plate 360, as shown in FIGURE 16C, of approximately the same thickness, 
incident angle, and glass type. However, correction plate 360 must be rotated 
90 degrees about axis Z with respect to dichroic filter 361- Correction plate 360, and 
dichroic filter 361 impart an equal but opposite amount of astigmatism in the 
10 transmitted wavefront, thereby canceling each other. Therefore, light striking 
detector 342 is free of astigmatism. This configuration leaves a small amount of 
^ residual coma. Yet, the optical performance is very close to the diffraction limit. Note 

%j that because light striking a detector 341 is reflected light, as opposed to transmitted 

CP light, there is less concern about astigmatism. 

; \ 15 Those of ordinary skill in the art will appreciate that the correction plate can be 

SI placed in many alternative positions, with adjustments in its thickness, material, and/or 

; s angle, relative to the propagation of the light. Note that FIGURE 16C shows only part 

;»i of an optical system, including only some of the elements employed after a collection 

u lens 340. Of course, the flow cells, optional light sources and additional detectors (see 

i 1* 20 FIGURE 1 6A) are included in a complete system. 

Any of the non-distorting spectral dispersing embodiments can be constructed 
using an additional objective lens 48 and slit 52, to form a confocal stop arrangement as 
shown in FIGURE 16B. FIGURE 16B also illustrates that detectors may be placed in 
either the transmission or reflection paths of dichroic filters. Either of the multiple 
25 detector embodiments may constructed such that the detectors receive light transmitted 
through the dichroic filters, reflected by the dichroic filters or in a combination of 
transmission and reflection as illustrated in FIGURES 16B. 

Note that FIGURE 16B also indicates optional ancillary components that might 
be beneficially incorporated into imaging systems in accord with the present invention, 
30 such as a pumping system, a reservoir for fluid into which encoded beads are entrained, 
and a reservoir for beads that have already been imaged and read. It should be 
understood that such components could be similarly incorporated in the other imaging 
systems of the present invention. 

In all spectral decomposition and imaging embodiments, it is anticipated that 
35 the magnification is set such that one detector pixel in object space is roughly the size 
of a single reporter. Given a detector pixel size of 10 microns and a reporter size of one 
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micron, the magnification would be configured to approximately 10X. likewise, if the 
reporter were 0.5 microns, the magnification would be set to approximately 20X. 

As disclosed in the above-referenced U.S. Patent No. 6,21 1,955 and as shown 
in FIGURES 1 and 4, multiple legs of the spectral decomposition and imaging system 
may be used to collect images and corresponding signals from multiple perspectives. 
Two legs are shown in FIGURE 1; however, if epi illumination is applied as shown in 
FIGURES 3 and 4, four legs may be used to view the beads from each side of the 
cuvette. Light emitted from the beads may also be collected with objective lenses that 
can be optically coupled directly to the cuvette to improve the numeric aperture. Jn 
addition, each optical leg can image beads at multiple focal planes to improve focus on 
reporter beads which may be out of focus at a different focal plane. 
Pixelated Detection 

Pixelated detection employs frame-based CCD image collection, in which a 
CCD detector views beads in flow, in a freeze frame fashion. This method requires the 
integration time to be very short to prevent blurring. A short integration time is 
achieved either with a strobed light source, or a continuous light source combined with 
a shuttered or gated detector. In either case, the short integration time reduces the 
signal-to-noise ratio and the ultimate sensitivity of the approach with fluorescence 
signals. Further, frame-based cameras require time to transfer data out of the pixelated 
detector, during which no images are acquired, and beads of interest can escape 
detection. However, these types of detectors are readily available, inexpensive and do 
not require an accurate knowledge of the velocity of the beads in flow. 
Alternative Approach for Pixelated Detection 

Another approach useful in bead imaging for pixelated detection also employs 
frame-based CCD image collection, but does not rely on strobed illumination or 
shuttered detection to freeze image motion. Instead, a rotating or oscillating mirror (not 
shown) is used to compensate for bead motion to produce a still image on the detector. 
This technique may employ continuous illumination, thereby achieving higher levels of 
sensitivity than strobed systems when analyzing fluorescence. However, this approach 
requires both an accurate measurement of bead velocity and a very stable fluid 
pumping system, since the inertia of the minor prevents compensation for rapid 
changes in bead velocity. 
A Further T echnique for Pixelated Detection 

Yet another method for pixelated detection uses TDI CCD image collection. In 
TDI detection, the electronic signal produced within the detector by an incident image 
is moved down the detector in synchrony with the motion of the image. In this manner, 
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signal integration times can be increased over conventional frame imaging modes by a 
factor exceeding 1000 fold. 
Bead Illumination 

As previously discussed in conjunction with several embodiments of die 
5 present invention and as shown in FIGURES 3 and 4, several different illumination 
systems may be employed to illuminate the beads in flow, A standard approach 
involves illuminating the beads in flow with a light source providing light along a path 
oriented orthogonal to the spectral decomposition and imaging system. Alternative 
modes of illumination, such as those shown and FIGURES 3 and 4 and disclosed in 
10 above-referenced U.S. Patent No. 6,21 1,955, allow for the generation of bright field, 
J3 dark field, phase contrast, fluorescence and EPI fluorescence imagery. U.S. patent 

application Serial No. 09/689,172, entitled 'Multipass Cavity for Illumination and 
Excitation of Moving Objects," filed on October 11,2000, discloses a method for 
illumination in which the number of photons incident on the beads may be increased by 
15 a factor of 1 0 or more. 
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\ j The design of the illumination system also allows the use of pulsed light source 

* or other strobed sources for high sensitivity fluorescence measurement without any 

I«! need to strobe in synchrony with bead flow. A high aspect ratio also allows for highly 

M efficient coupling of linear array diode illumination into the cuvette. 

™ 20 Bead Velocity Measurement 

For the second and third pixelated detection techniques, wherein an accurate 
knowledge of the bead velocity is required, either a frequency domain velocity 
measurement (FDVM) technique or a time domain velocity measurement (TDVM) 
technique, both as disclosed in US. patent application No* 09/939,292, filed on 
25 August 24, 2001 and entitled "Measuring the Velocity of Small Moving Objects Such 
as Cells" can be employed, In FDVM, a large FOV is imaged onto a ruling of opaque 
and transparent bars. Motion of the objects within the FOV causes modulation of their 
intensity as they pass across the ruling. The modulation frequency is proportional to the 
velocity of the objects and can be determined using Fast Fourier Transform analysis. 
30 Velocity can also be determined using two detectors in a conventional 

time-of-flight measurement scheme, though with very restricted throughput. However, 
time-of-flight measurements become more complex when throughput increases and the 
times of flight of multiple objects are measured simultaneously. Such systems can fail 
when correlation is lost between the entry and exit times of the objects in view. The 
35 time-of-flight system preferably used relies on an improved scheme wherein the 
waveforms produced by the entry and exit detectors are cross-correlated to detect phase 
changes that are indicative of changes in velocity. 


BlOLQ0n-H26/0UnAP.doc 


-32- 


Sample Handling 

A standard pumping technique of the type used in prior ait flow cytometers can 
be used to hydrodynamically focus the beads with micron-scale positional accuracy 
within the cuvette. Hydrodynamic focussing of the beads ensures they are located in or 
near the plane of focus of the optical system. Handling reporter labeled beads in 
suspension in the present invention eliminates the need to array beads in a monolayer 
on a substrate. When analyzing large samples containing a billion beads, the substrate 
methodology becomes impractical due to the number of substrates or the size of a 
single substrate required to carry a billion or more beads. A sample of four billion, 
10 micron beads occupies a volume of as little as 4.0 cm 3 or 40 milliliters. By 
comparison, when arrayed in a perfect tightly packed monolayer, that same four billion 
bead sample would occupy a square area of over 24 inches on a side. Or, when 
produced in more manageable format such as a standard microscope slide, the four 
billion bead sample will occupy more than 1728 slides. Other advantages of handling 
beads in suspension include the ability to sort individual beads at high speed, the ability 
to collect bead images from multiple perspectives, and the ability to hydrodynamically 
focus beads into a tight focussable core near the focal plane of the collection system. 
To increase throughput in flow, the sample stream may have a high aspect ratio in the 
plane perpendicular to flow. This technique is discussed in US. Patent No. 5,422,712 
"Apparatus for Measuring Fluorescent Spectra of Particles in a Row" the specification 
and figures of which are specifically incorporated herein by reference for purposes of 
providing background information regarding the technique commonly referred to as 
broad flat flow. Broad flat flow allows objects to flow in a plane as opposed to single 
file through the flow cell. This increases throughput by allowing beads to flow more or 
less parallel to each other while maintaining a tight plane of focus. Those skilled in the 
art will appreciate that broad flat flow can easily be created using commercially 
available flow cells as shown in FIGURE 16B containing flow cell 290, with a cross 
section that is elongated in an axis perpendicular to both the flow and optical axes. 
Spatial and Spectral Corrections to Pixelated Imagery 

In the second through fifth embodiments, a spatial and spectral correction may 
optionally be applied to the signals coming off the pixelated detectors) to improve the 
integrity of the decoding process. A method and apparatus for such a system is 
disclosed in US provisional patent application Serial No. 60/286,713 entitled 'Method 
and Apparatus for Cross-Talk and Spatial Registration Correction for Multi-channel 
Imaging." In this scheme, spatial registration errors between channels on the detector, 
or spatial registration errors between detectors, is corrected to more accurately 
determine the origin of light projected upon the detector. Likewise, once a spatial 
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correction has been applied, a spectral correction is applied to remove any spectral 
cross-talk between channels or detectors to more accurately determine the spectral 
content of light projected onto each pixel on the detectors). 

All figures illustrating imaging systems in the present invention employ 
pixelated detectors. In some embodiments multiple images of reporter labeled carrier 
beads are projected upon a single detector as illustrated in FIGURE 11. In other 
embodiments such as that illustrated FIGURE 16A multiple detectors are employed, 
each containing a different imagery of the same reporter labeled carrier beads. 

Ideally, each of the reporters or carriers being imaged by the present invention 
would deliver light comprised of wavelengths entirely bounded by the edges of the 
passband of a particular channel. In that case, each source type (bead fluorochrome) 
would appear in only one channel of a detector or on only one detector in the case of a 
multiple detector embodiment However, in many cases the fluorochromes used in the 
reporters, carriers or binding signals span a range of wavelengths broader than the 
passband of the associated imaging filter. This is clearly illustrated in FIGURE 23, 
where emission spectra of several fluorochromes are superimposed over the passbands 
of several channels on a wavelength versus intensity plot In this case, light from each 
source 201 will be received by two or three channels 202. The signal conveyed by each 
channel, then, is a composite of information from multiple sources. This can make 
reporter identification difficult, especially when reporters are intensity coded as well as 
spectrally coded. Therefore, the signals from the various detectors or channels must be 
operated on in a way to remove the crosstalk to accurately determine the spectral 
content and intensity for each reporter, carrier or binding signal 

Fundamentally, to remove crosstalk a set of linear equations as shown below 
must be solved . 

s x = + (Xy^ + a l3 m 3 

and 

where : 

rrtj = measurement from channel j , 
Sj = a characteristic parameter of source i, and 
= weighting coefficient, source i into channel j . 

The weighting coefficients carry information about the sources and about the 
channels used to collect the measurements. The equations are solved using the methods 
of linear algebra, wherein the variables % and mj are treated as vectors and variable 
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is treated as a matrix. The general tenia for this process is "crosstalk correction", since 
it removes the infomiation from a source that has fallen into a channel adjacent of the 
correct channel. This information spillover is analogous to crosstalk in bundled 
telephone lines, which results in a listener hearing the conversation from another line. 

This set of equations must be solved for each pixel location, ft is essential, 
then, that the images from all of the channels be precisely aligned with one another so 
that the correct measurements are entered into the equations for each pixel location. 
Therefore, a computational process is applied to the images belonging to a field of view 
before the crosstalk correction is applied. The "field of view" is the scene captured in 
an image. The images belonging to a field of view may be referred to as the "image 
ensemble" for that field of view. 

FIGURE 24 is a block diagram of the operations required to cany out the 
spatial and spectral corrections. It includes two significant image processing 
operations, applying spatial corrections executed in a block 305 and applying spectral 
cross talk corrections, executed in a block 306. In block 305, the images from a 
plurality of channels are spatially aligning, and in block 306, corrections are applied to 
remove the channel-to-channel crosstalk. There are two classes of information that 
support the operations executed in blocks 305 and 306. The first is the class of 
calibration constants, derived from on-line calibration of the instrument during its 
operation. The second class of supporting information is that for which stored tables of 
constants are accessed during operation, but not modified. 

In a block 301, a calibration image that is generated from a plurality of offset 
stripes (the spatial offset for the calibration image being known) is introduced to the 
instrument The instrument then computes the spatial offsets between the stripes in a 
block 302. The X,Y (horizontal and vertical) spatial offsets are determined in a 
block 303. Note that the image misalignment responsible for such spatial offsets may be 
subject to thermal drift and other factors that may vary during operation. The X,Y 
spatial offsets are used in block 305, along with other information, to apply spatial 
corrections. 

The stored constants will have been derived from measurements and from the 
known characteristics of the objects to be imaged. Different stored constraints will be 
applied to the spatial corrections executed in block 305 than are applied to the crosstalk 
corrections executed in block 306. For example, the general positions of the channels 
relative to each other known, and these positions are required for the spatial corrections 
executed in block 305. The stored data relating to the general positions of the channels 
relative to each other are provided to a processor executing the spatial corrections of 
block 305 in a block 308. With respect to the crosstalk corrections executed in 
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block306, stored data relating to the inverted source coefficients are provided in a 
block 309. It should be noted that additional stored measurements and known 
characteristics can also be incorporated into the second class, and provided to the spatial 
correction operation or the crosstalk correction operation, as appropriate. 

The image processing stages can be applied in real-time as images are collected 
during system operation, or by offline access to stored image files. The images to be 
spatially and spectrally corrected in blocks 305 and 306 are provided by a block 304. 
The resulting corrected images are output in a block 307, and are suitable for further 
processing and analysis, free of the information degradation caused by crosstalk and 
misalignment 

A block diagram of the general operations of the preferred embodiment of the 
present invention is shown in FIGURE 25. An "Image Ensemble" refers to a set of 
images, all of which depict the same field of view, but each of which has been 
constructed of signals from a particular channel. The number of images in an ensemble 
is equal to the number of channels in the instrument The image ensemble is generated 
in a block 401. 

The X (horizontal) and Y (vertical) offsets are required in order for the 
alignment process to operate on the image ensembles. These offsets are provided in a 
block 403. As discussed above with respect to FIGURE 24, calibration images are 
processed to determine these offset values. The calibration process may be ran as a 
preliminary step in the operation of the instrument The calibration may be repeated 
periodically to track drift in the image registration caused, for example, by changes in 
temperature. The process of generating the offsets is illustrated by a block diagram in 
FIGURE 26. 

In a block 404, images are aligned by shifting rows by pixel increments. Jn a 
block 406, an aligned image is generated by a floating point interpolation in a, using 
interpolation kernels generated in a block 405, based on the X.Y offsets determined in 
block 403. In a block 407, the crosstalk correction process described above is applied, 
using a spectral coefficient matrix provided by a block 408. Then in a block 409 the 
8-bit grayscale is restored, resulting in a final image ready to output in a block 402. 

Referring now to the block diagram of FIGURE 26, for each image, represent 
by a block 501, a first operation in the generation of the spatial offsets is that of 
boundary detection, executed in blocks 502 and 503. From the boundary information 
determined in block 502, a Fourier Transform is applied in a block 504, and then a Fast 
Fourier Transform is applied in the reference channel. From the boundary information 
determined in block 503, a Fourier Transform is applied in a block 505, and then a Fast 
Fourier Transform is applied in the data channel. The results of both Fast Fourier 
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Transforms are multiplies in a block 508, and an Inverse Fourier Transform is applied 
to the result in a block 509. Then accumulated 2D correlograms (1 per data channel) 
are generated in a block 510. The process is repeated for additional images. A detailed 
description of the process illustrated in the block diagram of FIGURE 26 is provided 
below. Note mat FIGURE 27 illustrates exemplary images before and after correction. 

A two-dimensional gradient operator is used to suppress flat surfaces and to 
enhance object boundaries. The operator builds signal in regions where the slope of the 
intensity is large in both the vertical and the horizontal directions. The linear 
expression for this operator is as follows: 

G= ±& 

dy dx. ' 

Object boundaries carry all of the energy in the images transformed by the 
gradient operator. What is needed next is a way to overlay each data image onto the 
reference image and measure how much vertical shift and horizontal shift is required to 
line up the object boundaries. Cross-correlation is a preferred method for measuring 
these shifts, or offsets, between images. 

The cross-correlation of two signals is comprised of convolving the two signals 
and extracting information from the output of the convolution operation, which is called 
the correlogram. The operation of convolution in the time domain is defined by the 
following equation: 

fi(t)®f 2 (t)=£f,(^)f 2 (t-X)dX. 

The value of the convolution for every time sample, t, is the sum over infinity of 
the product of the two functions, but with the second function offset by time L The time 
delay between the two signals can be determined by finding the peak in the correlogram. 

For image realignment, the convolution is applied in two dimensions to 
functions in x,y space, as follows: 

oo 

f t (x,y) ® f 2 (x,y) = JJf,(8,ll)f 2 (x -e, y-^dedi] . 

The alignment method is based on the premise that there will be similar 
structures in the images to be aligned. In the ideal case, the second image is identical to 
the first except for a shift in x and y. Furthermore, the shift can be represented as the 
convolution of the function with the twoKiirnensional Dtrac delta function, thus: 

Mx>y) = f 1 (x-x 0 ,y-y 0 ), 
and 

fii*,y) = /,(*, y)®S{x- y- y 0 ). 
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The image offsets, Xo and yo, can be measured from the shift in the two- 
dimensional correlogram, since: 

C u =f 1 (x,y)®f 2 (x,y) = f l (x,y)®f 1 (x,y)<8)6(x-x 0 ,y»y 0 )^ 

The Fourier Transfoim provides a convenient method for accomplishing the 
convolution of two images, as seen in the theorem: 

/iU y) ® f 2 (x, y) = F- x \F x {m x% & y ) • F 2 (o> x> a> y )] 
where : 

F(® x9 w y ) = Fourier Transform of f(x,y), 
F~ l [X(a> x ,a> y )] = Inverse Fourier Transform of X(a) x ,co y ). 

In the frequency domain, then, the two-dimensional convolution becomes a 
matrix multiplication of the spectra of the two images. The two-dimensional 
correlogram is obtained by performing the inverse Fourier Transform of the product. 

Because objects may occupy only a few pixels, and because the resolution of 
the imaging system may be on the order of a pixel width, alignment of one image to the 
next to a resolution of a fraction of a pixel width is necessary for accurate crosstalk 
correction. The true peak of the correlogram will rarely be centered on a pixel. The 
true peak can be located, however, by analyzing the shape of the region around the 
pixel of maximum amplitude. If an accurate equation can be defined for the 
correlogram amplitude in that region, the true location of the peak can be found to a 
resolution better than a single pixel width. 

For cross talk correction an image is reconstructed for alignment by application 
of two-dimensional interpolation. In this process, once the image has been aligned with 
the reference image to the nearest pixel, the new amplitude value for each pixel is 
computed as the weighted sum of a group of surrounding pixels. The values of the 
weighting coefficients for the interpolation depend on the amount of shifting to be 
accomplished in the vertical direction and the horizontal direction. These shift 
distances will be less than one pixel width. 

With aligned data images in place, the crosstalk correction can be 
accomplished. FIGURE 27 illustrates imagery before and after cross talk correction. 
Imagery 1601 clearly shows objects imaged upon a detector in the presence of 
uncorrected cross talk. Imagery 1602 illustrates the same imagery after cross talk has 
been removed Depending upon the means for used for tagging reporters, 
fluorochromes, quantum dots etc., it may or may not be necessary to remove cross talk 
from imagery projected upon the detectors) before decoding. 
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Analvzing Reporter labeled Beads Using the Flow Imaging System of the Present 
Invention 

Decoding of reporter labeled beads by the second embodiment of the flow 
imaging system is depicted in FIGURE 17. Each bead and its associated reporters are 
5 simultaneously imaged multiple times across several rows covering a horizontal region 
on the detector, each image corresponding to a different emission, reflectance or 
absorbance band. In FIGURE 17, the solid substrate, or carrier bead 346, is indicated 
by the large circle, while reporters 348a, 348b, 348c and 348d are depicted as a smaller 
circle bound to the larger bead, 

10 Each of these reporters can be seen in the (Optional) reference channel of the 

detector. The imagery in this section can be generated by light source scatter signals 
from the beads and reporters or by absorbance signals generated from band limited 
brightfield illumination. For reporters using fluorescence to uniquely identify their 
optical signature, the reference channel is optional if each reporter contains at least one 

15 type of fluorochrome. A binding signal is indicated by a filled circle 346b in the 
binding signal section of the detector. A signal in this channel indicates that an analyte 
has bound to the bead as aresult of the assay. Further, the intensity of the signal in the 
binding signal channel can be used to further characterize the binding event to 
determine the quantity of the analyte present in the assay. 

20 Reporter signals are imaged in the remaining blue, green, yellow and red 

channels of the detector. The present invention maintains the relative positions of all 
reporter imagery in each channel Thus reporter 348a is associated with imagery 342a 
in the yellow channel and imagery 342b in the red channel of the detector. 
Reporter 348a contains both yellow and red fluorochromes, and since no other imagery 

25 appears in the blue or green channels, no other fluorchromes are present in 
reporter 348a Similarly, reporter 348b is associated with imagery 343a in the blue 
channel and imagery 343b in the red channel, reporter 348c is associated with 
imagery 344a in the green channel and imagery 344b in the yellow channel, and 
reporter 348d is associated with imagery 345a in the blue channel and imagery 345b in 

30 the green channel 

Such reporter imagery must be matched to a known reporter legend, to identify 
an object that is associated with a specific reporter. It should be noted that reporter 
labeled beads can be associated with many different types of objects, such as cells, 
chemical or biological compounds. As discussed in the Background of the Invention, 

35 certain types of DNA related research involves the preparation of libraries of 
synthesized compounds, such as an oligo library generated using the SAP technique 
described above. Such compounds are produced by combining a plurality of subunits 
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(for example, nucleotide bases) together to form a larger compound (such as an 
oligomer). While several different methods can be employed to produce such 
compounds, if each subunit is individually tagged with an identifying characteristic, 
such as a fluorchrome, then the presence of a reporter encodes the identity of the 
subunit in the compound, and its location relative to other subunits, based on a known 
reporter legend FIGURE 19 illustrates a reporter legend, that was generated when a 
oligo library comprised of nucleotide bases was generated. According to the reporter 
legend of FIGURE 19, the presence of the yellow-red reporter 348a indicates 
cytosine(Q in a first position in the oligo bound to bead 346. Hie presence of the 
blue-green reporter 348d indicates an adenine (A) in the last position of the oligo bound 
to bead 346. In this manner, all positions of an oligo can be determined by identifying 
the color combinations present in each reporter. Using the present invention it is very 
straight forward to decode the XGTA" oligo sequence synthesized on bead 346. 
Furthermore, the presence of the binding signal (see FIGURE 17) indicates whether or 
not the complimentary "GCAT" sequence has bound to the bead during the assay. It 
should be noted that while reporter labeled beads will be useful in identifying all the 
subunits of a compound that include individually labeled subunits, it is anticipated that 
it will also be useful to label (and read) compounds not at the subunit level, but just at 
the compound level. Thus the present apparatus and method is not limited to imaging 
and reading only compounds that include individually labeled subunits. 

The decoding process described above is identical for the fourth and fifth 
embodiment of the flow imaging system described above, but each spectral band is 
imaged onto a separate detector. Although this embodiment requires more detectors, 
each detector can have relatively few pixel columns since there is only one spectral 
band per detector. This embodiment is also more optically efficient because the light 
reaching the detector passes fee filter set once, unlike the previous embodiment where 
light must pass through each filter twice before exiting the dispersing component and 
reaching the detector. 

The decoding process used for the first embodiment is similar in concept to the 
process employed for the second through fourth embodiments, except that the reporter 
images must be generated by first deconvolving emission spectra of the fluorochromes 
used from the imagery on the detector. Deconvolution of the emission spectra fiom the 
bead image requires significant image processing, which is a disadvantage for real-time 
analysis. However, the advantage of this first embodiment is that it potentially provides 
the highest spectral discrimination power, but uses only one detector and one optical 
dispersing element, reducing system cost 
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FIGURES 18A and 18B include flowcharts 310 and 410 that provide the series 
of logical steps employed to read reporter labeled beads in accord with the present 
invention. Referring to FIGURE 18A, in a start block 312, a first bead is read. In a 
block 313, spectral and/or spatial corrections are optionally applied to the pixelated 
5 signals from the detector to more accurately determine the spectral content in each pixel 
of the imagery. Next, the data from the bead imaging are segmented into a plurality of 
different reporter signature channels in a block 314. The number of signature channels 
in the present invention is chosen to be equal or greater than the number of different 
colors present in any reporter. In FIGURE 19, the reporter legend is based on reporters 
1 0 that comprise no more than 4 colors. Thus, FIGURE 17 includes four separate reporter 
signature channels, including "BLUE," "GREEN" "YELLOW " and "RED" 

In a block 3 16, the deconvolution step discussed above is performed, but only if 
required, as would be the case for the first embodiment of the imaging and spectral 
fif I decomposition. In a block 3 1 8 the binding signal channel is checked to determine if an 

| 1 5 analyte has bound to the bead indicating the presence of the analyte in the sample being 
£i tested If no binding signal is present the logic determines if more beads are to be read 

in a decision block 324. If in decision block 324 it is determined that no more beads are 
to be read, a list of bead identities and binding strengths are generated in a block 326, as 
;!l will be described in more detail below, with respect to FIGURE 18B. If in decision 

ry 20 block324 it is determined that more beads are to be read, the next bead is read in a 
;« : 5 block 328, and the logic loops back to block 3 13 as described above 

If in decision block 318, it is determined that a binding signal is present, the 
presence and or strength of the binding signal is stored in a memory in a block 320. In a 
block 322, the bead is decoded as described in detail below in conjunction with 
25 FIGURE 18B. In the decoding process, when a valid bead signature is found, the 
signature or compound identity is stored along with the binding indicator. At this point 
the logic determines if more beads are to be read, in decision block 324, as described 
above. 

The method for determining compound identity is illustrated by flowchart 410 
30 in FIGURE 1 8B, which as noted above is executed in block 322 of FIGURE 18A, as 
noted in a block 420 in FIGURE 18A. After bead imagery has been segmented in each 
channel, individual reporters are evaluated to determine their identity. In a block 422 
the imagery corresponding to each reporter is located in each signature channel. This 
may be accomplished by any means, including but no limited to, evaluating imagery in 
35 each signature channel to find the top left most reporter image and its corresponding 
imagery in other channels. Moving from left to right and top to bottom across the 
segmented region in each channel, every reporter is located in each channel. In a 
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block 424 a reporter is evaluated for its spectral content and the results are stored In a 
decision block 428, the logic determines if more reporters are to be analyzed. If so, the 
logic loops back to block 424 to repeat the process for each reporter found in the image. 
Once it is determined in decision block 426 that no more reporters are to be analyzed, 
then in a block 428 the reporter signatures are compared against the reporter legend of 
FIGURE 19 to determine the compound identity. 

In a decision block 430, the logic determines if the bead has a valid signature. 
If the signature is not valid, that occurrence is noted in a block 432. If the signature is 
valid, the signature is stored in a block 434. A signature is invalid if it is not supported 
by the reporter legend -The process illustrated in FIGURE 18B is repeated for each 
bead. It should be understood that when such a process is automated and executed by a 
computing device, the process of analyzing a single bead occurs very rapidly. 

Flow Analysis of libraries for DN A Sequencing, Polymorphism Analysis, and 
Expression Analysis 

Given a bead library that comprises a complete set of DNA oligomers of length 
N, a subset of the library will hybridize to any piece of single stranded target DNA of 
approximately length N or greater. If the complementarity of the target DNA to a bead 
is less than iV bases, hybridization strength (measured by the melting temperature, 7*J 
will drop. Six reporter colors employed in a binary encoding scheme are sufficient to 
encode a bead library of oligos up to length sixteen. As discussed above, the choice of 
oligo length affects the number of beads necessary for a complete library, the number 
of reporters necessary to encode the library, and the strength of hybridizatioa Shorter 
oligos require the smallest bead sets, the fewest reporters, and have the lowest tolerance 
for hybridization mismatch. 

The shortest practical oligo length is approximately a 10-mer. A complete 
10-mer library requires approximately one million beads, which can be analyzed in 
under an hour at a rate of a few hundred beads per second The T m of a perfect 10-mer 
complement is near room temperature, particularly for ATI rich sequences, so just a few 
base mismatches will likely cause the target DNA to dissociate from a bead at room 
temperature. Oligos of length 1 1 through 16 have melting temperatures significantly 
higher than room temperature and can therefore tolerate more mismatches than oligos 
of length 10, given the same reaction stringency. For the purposes of the following 
description, it will be assumed that a bead library of all possible 10-mers is used. 

If a sample of target DNA consisting of multiple identical copies of length M is 
hybridized to a complete library, a subset of the beads numbering approximately MAO 
will hybridize at random to the different strands of target DNA based on the sequences 
of their attached oligos. Each of the hybridizing oligos will have sequences that are 
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offset by one nucleotide, forming a contiguous overlapping string of sequences, as 
shown in FIGURE 20. Labeling the target DNA directly or indirectly with at least one 
luminescent probe will cause the bead to which the DNA is bound to emit light with the 
characteristic emission spectrum of the probe* If the probe spectrum is distinct from the 
5 spectra used for the reporters, the hybridization can be detected as shown in 
FIGURE 21. 

After hybridization, the bead library is analyzed by flow imaging and the beads 
that emit a binding signal are decoded to reveal their associated oligonucleotide 
sequences. These sequences are then assembled into a "contig?' of overlapping 
10 sequence by starting with a single oligo sequence and searching the remainder of the set 
^ of oligo sequences that bound to the target to find the closest match. The search can be 

!n highly constrained by the knowledge that all but a few of the bases near the end of the 

SI next oligo to be added to the contig must overlap the end sequence of the contig. The 

j sequencing capacity of a bead library is constrained by the probability that the contig 

i jj 15 will branch and loop during the sequencing of a single stretch of DNA or that it will 
Si cross multiple fragments that are being sequenced simultaneously. The maximum 

f , length of a DNA sequence is approximately the square root of the size of the bead 

q library. Hence, a complete 16-mer library of over four billion beads can sequence a 

H stretch of DNA approximately 65,000 bases long. 

^: 20 As illustrated in FIGURE 20, the single-base offset of each oligo in the contig 

U assures that there is AT-fbld sequence redundancy everywhere but at the ends of the 

contig. Further, the bead library can embody redundancy by including multiple copies 
of each unique bead (and corresponding oligo). Having multiple beads of each type 
enables the system to tolerate some fraction of unanalyzed beads due to system latency, 
25 defocus, out-of-view reporters, etc. 

The sequencing protocol described herein can be implemented with a single 
target DNA sample or a pool of multiple targets, which will result in multiple contigs. 
The target DNA can be from any source, including genomic DNA, viral DNA, DNA 
produced by an amplification reaction such as a Polymerase Chain Reaction (PCR), or 
30 any other source known to those of ordinary skill in the ait In addition, DNA produced 
by an RNA template, RNA itself, or any other molecule known to hybridize to DNA, 
such as a peptide nucleic acid, can be used. 

FIGURE 21 illustrates the application of the flow imaging sequencing method 
to determine the identity of SNPs, generally as indicated in the steps of FIGURE 22. A 
35 target DNA molecule is shown with a polymorphic site that can contain either a G or a 
T component. In sequencing the target, the ten beads that hybridize across the SNP site 
will either have a C or an A component at the SNP position if the source DNA is 
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homozygous. If the target DNA is heterozygous, the beads binding across the SNP site 
will consist of two populations with either a C or an A component at the SNP position. 
If the DNA used as a template for the amplification of the SNP regions is pooled from 
numerous sources, the relative frequency in the population of the different 
polymorphisms at a given SNP position will be reflected in the relative abundance of 
the amplification products. The relative polymorphism frequency will therefore also be 
reflected in the binding signal intensity ratios of the different bead sequences that span 
the SNP site. Numerous polymorphic sites in a large genome can be analyzed 
simultaneously by amplifying and sequencing short regions from the genome that 
incorporate the polymorphisms. For example, the 50 nucleotides flanking an SNP 
require approximately 40-50 distinct beads for analysis. Since this approach is 
equivalent to DNA sequencing, the library capacity constraints are similar. A million 
bead library of 10-mers can therefore analyze the relative frequency of a few tens of 50 
base SNP fragments. If 12-mers are employed, the library size grows to 16 million, 
enabling the detailed analysis of roughly 100 SNPs. 

An alternative means for analyzing SNPs is to employ polymorphically 
selective amplification during the production of the target DNA. That is, if a specific 
polymorphism is present, it will be amplified from the SNP region. Absence of the 
polymorphism results in no amplification. Because the polymorphic selectivity is 
inherent in the amplification process, the bead analysis serves only to identify the 
presence or absence of a reaction product. The presence of a reaction product can be 
determined by a binding signal on as few as one unique bead associated with the 
product, which is a far more relaxed constraint on the bead library's analytical capacity 
than the case of DNA sequencing and greatly increases the analytical capacity. For 
example, if 50 base stretches of DNA are amplified from multiple SNP sites and pooled 
for analysis, each SNP fragment must differ from all the others by at least one bead. 
Hence, there can be extensive cross-branching of the SNP fragment hybridization 
patterns. The analytical capacity of a library therefore approaches the size of the library 
divided by the number of beads needed to hybridize to a given SNP fragment In the 
case of a million bead 10-mer library with 50 base fragments, on the order of 10,000 
SNPs can be analyzed simultaneously, 

A variation on the above method of SNP genotyping is to employ a different 
probe color for each base that may be present at an SNP site. In this manner, relative 
polymorphism abundance can be obtained without the need to sequence the SNP 
regions. 

In addition to DNA sequencing and polymorphism analysis, flow imaging of 
bead/oJigo libraries can be used to analyze gene expression. In the case of gene 
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expression analysis, specific sequences from each of the genes of interest are extracted 
from RNA, Extraction can take the form of splicing out and labeling sections of RNA 
itself or converting the RNA to labeled DNA, for example by reverse transcriptase 
PCR. The same basic principle of the second method of polymorphism analysis applies 
5 to expression analysis, namely that there is no need to sequence the entire expressed 
gene. Instead, the presence of a specific fragment, as indicated by the binding of at 
least one unique bead to the amplified expression sequence, indicates the expression of 
a particular gene. As in sequencing polymorphic sites, the intensity of the binding 
signal can indicate the level of expression. 

10 The hybridization of DNA to bead-bound oligos is reversible, raising the 

possibility that a bead library could be used a number of times. There are several ways 
to restrict the reuse of a bead library based on the fact that the correspondence between 
each bead signature and its associated oligonucleotide sequence is determined and 
controlled during the production of a bead lot and can be varied from lot to lot The 

15 information necessary to decode oligo sequences from the bead signatures detected 
during analysis can be restricted to prevent the unauthorized use of a bead library. The 
analysis platform could identify the unique lot number of the bead library via a barcode 
or other identifier and employ a corresponding decoding file distributed by the 
authorizing body to perform the analysis once, or a given number of times. Another 

20 method of restricting authorization might require the user's analytical platform to contact 
the authorizing entity with the bead lot number at the time of analysis, via the Internet 
for example, whereupon the identity and validity of the user is determined. Alternately, 
the beads could be physically altered during analysis, for example by illumination with 
260 nm light, to destroy the bound oligonucleotides and prevent reuse. 

25 The use of reporter labeled libraries in flow imaging is not limited to DNA and 

RNA analysis. Any molecule that can be synthesized on a bead, such as amino acids or 
drug candidates, can be encoded and read in flow. The concept can even be extended 
to living cells. Cells of different types can be encoded by attaching labels to their 
membranes or locating those labels within the cell itself. For example, a reporter with a 

30 unique signature identifying a corresponding specific receptor-specific antibody could 
target any cells expressing that receptor on the cell surface. In some cases, subsequent 
endocytosis of the receptor-antibody complex would then internalize the reporter bead 
as well. This binding would identify those cells expressing a specific receptor. 
Reporters could be internalized by the cell via a variety of mechanisms, including, but 

35 not limited to, phagocytosis, endocytosis, lipophilic conjugation, protein transduction, 
or alternatively through pores created in the cell membrane via exposure to toxins, such 
as streptolysin 0,tetanolysin, E. coli hemolysin, via electroporation, or using other 
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means known to one skilled in the art. Specific reporters can be used to indicate cell 
type, added over time to indicate cell age, and or incorporated upon the cell's exposure 
to certain drugs or environmental conditions. For encoding even greater diversity of 
cell types or cell processes, the concept can be extended to tagging cells with reporter 
5 labeled carrier assemblies, In this case, the cell acts as a carrier bead which carries 
other reporter labeled carrier beads. Note also reporter labeled compounds are not 
limited to multi-component compounds being labeled as there are synthesized Single 
component and already synthesized multi component compounds can be bound to 
beads . Reporters can then be used to label any such compound associated with a bead. 
10 Referring once again to FIGURE 22, a flowchart 350 provides the series of 

Q logical steps employed to generate and use an encoded bead library for either 

~? genotyping, sequencing, or expression profiling. In a block 352, a library of N-mer 

oligos is generated The combinatory SAP scheme described above can be beneficially 
W employed as well as a directed synthesis. While most of the current discussion has 

I ]\ 15 been directed to optical reporters that use color as a distinguishing characteristic, it 

Sit 

\a should be noted that other reporters, with other distinguishing characteristics, as noted 

above, can also be beneficially employed to generate an encoded bead library of the 
desired size. In a decision block 354, a user must select the type of analysis to be 
performed (genotyping, sequencing, or expression profiling). 

II 20 If genotyping is selected the logic proceeds to a block 356, and the genomic 
t?l DNA is amplified using primers for polymorphic regions of interest, as discussed 

above. If in decision block 354, expression profiling was selected as the analysis of 

choice, then in a block 358, RNA is amplified using primers for genes of interest, also 

as generally discussed above. Alternatively, if sequencing is selected in decision 
25 block 354, then in a block 360, a user would amplify genomic DNA using primer for 

extended sequences of interest 

Regardless of the analysis selected the next step is to hybridize the amplified 

DNA/RNA to the bead library in a block 362. fa ablock364, the beads are analyzed 

preferably by employing one of the flow imaging systems discussed above, to identify 
30 oligo sequences of positive beads. Once the positive beads are imaged the final step is 

to construct sequence contigs to identify genomic DNA sequences, polymorphic 

alleles, or expressed genes, in a block 366. 

Although the present invention has been described in connection with the 

preferred form of practicing it and modifications thereto, those of ordinary skill in the 
35 art will understand that many additional modifications can be made to the present 

invention within the scope of the claims that follow. Accordingly, it is not intended that 
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the scope of the invention in any way be limited by the above description, but instead 

be determined entirely by reference to the claims that follow. 

Flow Analysis of Libraries for Pharmaceutical Chemical Compound Library Analysis 

Contemporary pharmaceutical research often involves the discovery of novel 
ligands which bind to specific cell receptors, thereby producing a pharmacologically 
beneficial response. The ligands might consist of small peptides or other oligomers of 
a size and molecular weight that are potentially useful as pharmaceutical products. In 
the present invention, a large library of compounds may be created and attached to 
reporter labeled beads. This may be accomplished simultaneously as the compounds 
are synthesized through an SAP process, directed synthesis or binding pie-synthesized 
compounds to pre-encoded reporter labeled bead assemblies. In the present invention 
the entire library can then be exposed to numerous sets of receptor targets in an assay in 
which the targets are fluorescendy labeled, undergo a competitive binding assay, or 
other assay techniques known to those of skill in the art. The exposure can occur in a 
controlled fashion in a micro-titer plate well or in a micto-fluidic device. Preferably, 
numerous copies of the compound library can be pipeted into a micro-titer plate well, 
along with the multiple copies of a target After sufficient exposure time, the contents 
of the well can be withdrawn and analyzed (or analyzed in situ), in accordance with the 
present invention, to determine which compounds demonstrate binding affinity to the 
target The process can be repeated for multiple wells containing different targets. 
Alternatively, if the targets are labeled with distinct fluorochrames, several targets can 
be interrogated simultaneously in a single well with the compound library. In this case, 
one of the channels shown in FIGURE 17 may serve a dual purpose, containing both 
small reporter imagery superimposed over a larger image of the carrier bead, thereby 
indicating the presence of a binding signal. 
Flow Analysis of Reporter Encoded Cells 

In a fashion similar to that disclosed above, the present invention can be used 
with cells that are exposed to different compounds, viruses, bacteria or environmental 
conditions. During or preceding exposure, the cells can be labeled with reporters to 
encode the exposure history, and or the cell type used in an assay. In accordance with 
the present invention, this process can be conducted in wells of a micro-titer plate, 
within a iracro-fluidic chip, or on a solid support such as a microscope slide. Once the 
exposure cycle has been complete the cells can be analyzed to characterize relevant 
biological affects. FIGURE 17 illustrates how these cells would appear when projected 
onto a detector. In this case elements 346 and 300 represent imagery of the cytoplasm 
and or nuclei. 
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Although the present invention has been described in connection with the 
preferred form of practicing it and modifications thereto, those of ordinary skill in the 
art will understand that many other modifications can be made thereto. Accordingly, it 
is not intended that the scope of the present invention in any way be limited by the 
above description. 
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